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Abstract 

The  effects  of  microcracking  on  the  mechanical  properties  of  Salem  lime¬ 
stone  were  investigated  in  three  phases:  introduction  of  quantifiable  levels 
of  microcracks  by  thermal  treating,  mechanical  testing  of  limestone  sam¬ 
ples  with  varying  levels  of  microcracks,  and  modification  of  a  numerical 
model  to  incorporate  the  measured  effects.  Computed  tomography  scan¬ 
ning,  scanning  electron  microscopy,  and  optical  microscopy  (OM)  were 
used  to  observe  microstructural  changes  caused  by  the  heat  treatments. 
Mechanical  testing  was  performed  to  characterize  the  mechanical  response 
of  the  intact  and  damaged  limestone.  Quasi-static  tests  included  uniaxial 
compression,  triaxial  compression,  hydrostatic  compression,  and  uniaxial 
strain  /  constant  volume  tests.  Microcracking  did  not  affect  the  limestone’s 
strength  at  pressures  greater  than  10  MPa.  Dynamic  tests  were  performed 
using  a  modified  split  Hopkinson  pressure  bar. 

The  results  of  the  mechanical  tests  were  used  to  modify  the  HJC  model. 
Modifications  were  made  to  account  for  shear  modulus  degradation  and 
failure  surface  changes.  The  original  and  modified  HJC  models  were  used 
in  a  numerical  analysis  of  the  mechanical  tests  performed  in  this  work. 

The  modified  HJC  provided  better  results  for  damaged  material  when 
compared  with  the  quasi-static  and  dynamic  experiments.  This  work 
demonstrated  that  this  approach  is  useful  for  examination  of  the  effects  of 
microcracking  on  quasi-brittle  materials  and  can  be  used  to  improve  the 
predictive  capabilities  of  material  models. 
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1  Introduction 

1.1  Background 

This  report  examines  three  topics:  the  quantification  of  induced  mi¬ 
crocracking  in  Salem  limestone,  the  effects  of  microcracks  on  its  mechani¬ 
cal  response,  and  the  numerical  modeling  of  that  mechanical  response. 

The  first  two  topics  deal  with  techniques  to  induce  and  quantify  various 
levels  of  microcracking  and  determining  changes  in  the  mechanical  re¬ 
sponse  resulting  from  the  microcracks.  These  topics  are  aimed  at  increas¬ 
ing  the  understanding  of  the  effects  of  microscale  properties  on 
macroscale  material  response.  The  last  topic  introduces  modifications  to  a 
commonly  used  numerical  model  intended  to  improve  its  predictive  capa¬ 
bilities  for  materials  containing  defects  such  as  microcracks.  These  studies 
demonstrate  the  importance  of  material  characterization  for  varying  levels 
of  defects  and  how  that  characterization  can  be  used  to  improve  current 
predictive  models. 

The  material  used  in  this  study  is  a  quasi-brittle  rock  called  Salem  lime¬ 
stone.  Salem  limestone  was  chosen  for  this  work  mainly  because  it  is  a  rel¬ 
atively  homogeneous  material  and  its  quasi-brittle  behavior  should  be 
analogous  to  other  quasi-brittle  materials  such  as  concrete.  Its  relative  ho¬ 
mogeneity  simplifies  the  microscale  quantification  and  mechanical  testing 
by  reducing  the  number  of  microscale  mechanisms  involved  in  the  mate¬ 
rial  characterization  and  providing  more  consistent  results  to  make  the  de¬ 
termination  of  microcrack  effects  more  straightforward.  Characterizing  the 
effects  of  microcracking  on  this  limestone  should  provide  an  easier  transi¬ 
tion  to  performing  similar  work  on  more  complicated  materials  such  as 
concrete,  which  have  more  complex  chemistry  and  microstructure  than 
limestone. 

Understanding  the  mechanical  behavior  of  quasi-brittle  materials  is  the 
foundation  of  many  civil  engineering  applications  due  to  their  wide  use  in 
civil  and  military  structures.  These  materials  include  natural  geological 
materials,  such  as  limestone  and  granite,  and  engineered  materials,  such 
as  concrete  and  ceramics.  The  evaluation  of  the  mechanical  properties  of 
these  materials  has  been  largely  based  on  their  behavior  in  pristine  condi¬ 
tion,  and  that  information  is  used  when  designing  structures  created  from 
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these  materials.  However,  the  mechanical  behavior  of  these  structures  un¬ 
dergoes  significant  changes  during  their  life  cycle  due  to  ordinary  wear,  ve¬ 
hicular  impact  (in  civil  applications),  blast  wave  impact,  and  projectile 
penetration  (in  military  applications)  among  other  things.  The  changes 
that  the  structures  undergo  during  these  events  should  be  taken  into  ac¬ 
count  if  their  continued  use  following  these  events  is  anticipated. 

Vehicular  and  blast  wave  impact,  projectile  penetration,  and  even  ordinary 
use  causes  microscopic  damage  in  brittle  materials,  such  as  rock  and  con¬ 
crete,  in  the  form  of  microcracks.  Furthermore,  microcracks  can  be  caused 
by  stresses  from  mechanical  or  tectonic  loading,  differential  expan¬ 
sion/contraction  between  grains  from  thermal  loading,  expansion  of  fluids 
in  void  space  during  freezing/thawing,  or  by  various  chemical  processes. 
The  nucleation  and  growth  of  microcracks  significantly  affect  the  mac¬ 
roscale  behavior  of  the  material  by  reducing  the  elastic  properties,  induc¬ 
ing  material  anisotropy,  and  causing  irreversible  strains  due  to  crack 
openings.  A  better  understanding  of  microstructural  changes  in  brittle  ma¬ 
terials  caused  by  these  events  and  its  effects  on  the  macroscale  mechanical 
response  is  needed  to  more  accurately  model  the  damaged  response  of  a 
brittle  material. 

If  a  correlation  between  the  microstructural  changes  in  a  material  and  its 
effects  on  the  macroscale  response  is  determined,  it  can  be  used  in  several 
ways: 

1.  Quantification  of  microscale  changes  following  a  loading  event,  such  as 
a  projectile  impact,  can  be  used  to  predict  the  material  state,  e.g.  spa¬ 
tially  varying  microcrack  densities,  for  similar  events. 

2.  The  correlation  between  the  microstructural  changes  and  macroscale 
response  can  be  used  in  phenomenological  models  based  on  microscale 
phenomena  to  predict  the  material  state  following  a  loading  event  in 
lieu  of  actual  experimental  data. 

3.  The  experimental  data  for  the  macroscale  response  of  materials  with 
varying  levels  of  microcracks  can  be  used  to  establish  basic  principles 
to  guide  the  implementation  of  damage  in  empirical  models. 

This  report  addresses  these  issues  by  quantifying  various  levels  of  induced 
microcracks,  determining  the  changes  in  mechanical  behavior  for  various 
stress  states  and  strain  rates  due  to  the  microcracks,  and  applying  the 
principles  learned  from  these  tests  to  a  material  model. 
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1.2  Organization 

This  report  covers  the  three  main  topics  in  five  chapters:  the  quantification 
of  induced  microcracking  in  limestone  (Chapter  3),  the  effects  of  mi¬ 
crocracks  on  its  mechanical  response  (Chapter  4),  and  the  numerical  mod¬ 
eling  of  that  mechanical  response  (Chapter  5).  Chapter  2  provides  the 
current  state  of  each  of  these  topics  and  Chapter  6  summarizes  the  previ¬ 
ous  chapters  and  provides  recommendations  for  future  work. 

In  Chapter  2,  the  current  state  of  microcrack  characterization  following 
mechanical  and  thermal  loading  will  be  reviewed.  This  review  will  concen¬ 
trate  on  experimental  work  related  to  microcrack  evolution  instead  of  mi¬ 
crocrack  evolution  theories.  Following  the  microcrack  evolution 
discussion,  two  types  of  phenomenological  material  models  currently  be¬ 
ing  used  to  predict  the  behavior  of  damaged  brittle  materials  will  be  dis¬ 
cussed.  This  discussion  will  provide  details  on  how  the  models  implement 
damage  and  its  effects  on  the  macroscale  mechanical  response. 

Chapter  3  will  provide  the  experimental  techniques  used  in  this  work  to 
quantify  microcracking.  The  chapter  is  divided  into  four  sections  to  pro¬ 
vide  information  on  the  material  used  in  this  work,  the  experimental 
method  used  to  induce  damage,  the  experimental  methods  used  to  quan¬ 
tify  the  induced  damage,  and  the  results  obtained  from  those  experimental 
methods.  The  sample  preparation  methods  for  each  experimental  method 
used  will  also  be  given. 

In  Chapter  4,  test  methods  to  get  the  mechanical  properties  of  the  material 
will  be  discussed  along  with  the  results  of  those  tests.  Mechanical  testing 
was  performed  to  characterize  the  mechanical  response  of  the  material 
with  varying  amounts  of  microcracking  at  low  and  high  strain  rates  for  var¬ 
ious  loading  states.  The  sample  preparation  methods  are  given  along  with 
the  servo-hydraulic  loading  techniques  used  in  quasi-static  tests  and  the 
split  Hopkinson  pressure  bar  technique  used  in  the  dynamic  tests.  Chapter 
5  presents  a  numerical  analysis  of  the  mechanical  tests  discussed  in  Chap¬ 
ter  4.  The  numerical  analysis  was  performed  using  a  commonly  used  con¬ 
crete  plasticity  model.  Modifications  to  the  model  are  made  to  account  for 
changes  in  the  mechanical  behavior  due  to  the  induced  microcracking.  The 
results  from  the  original  and  modified  model  are  compared  to  demonstrate 
the  improvements  due  to  the  modifications. 

Lastly,  Chapter  6  is  the  conclusion  to  this  report.  This  chapter  summarizes 
the  previous  chapters  and  provides  recommendations  for  future  work. 
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2  Current  State 

In  this  chapter,  an  overview  is  presented  of  the  current  state  of  the  three 
topics  examined  in  this  report:  the  quantification  of  induced  microcrack¬ 
ing,  the  effects  of  microcracking  on  the  mechanical  response  of  quasi-brit- 
tle  materials,  and  the  numerical  modeling  of  the  mechanical  response. 
Regarding  the  first  two  topics,  the  literature  most  relevant  to  this  work  will 
be  discussed  and  how  this  work  seeks  to  further  the  understanding  of 
these  topics.  For  the  numerical  modeling,  the  goal  is  not  to  give  an  exhaus¬ 
tive  list  of  the  available  models  for  quasi-brittle  materials  but  to  give  an 
overview  of  the  different  methods  to  simulate  damaged  mechanical  behav¬ 
ior  proposed  in  literature  along  with  an  in-depth  description  of  a  com¬ 
monly  used  material  model  that  is  modified  in  Chapter  5  based  on  the 
experimental  data  found  in  this  work. 

2.1  Background 

The  quasi-brittle  material  used  in  this  study  is  commonly  referred  to  as  In¬ 
diana,  Salem,  or  Bedford  limestone.  This  particular  rock  was  chosen  for 
this  work  because  of  its  relative  homogeneity,  the  availability  of  test  data 
from  other  work  (Cummins  1991;  Frew  2001;  Frew  et  al.  2001),  and  its 
mechanical  behavior  should  be  analogous  to  other  quasi-brittle  materials, 
e.g.,  concrete,  ceramics,  and  other  rocks.  These  materials  are  referred  to  as 
quasi-brittle  because  they  exhibit  moderate  strain  hardening  prior  to  at¬ 
taining  their  strength,  followed  by  a  strain  softening  behavior  (Karihaloo 
1995)- 

Generally,  quasi-brittle  materials  are  composed  of  a  hard,  solid  material 
(grains,  aggregates,  crystals)  and  a  softer  matrix  material,  or  bond  system. 
The  matrix  material,  even  in  intact  samples,  usually  contains  an  abun¬ 
dance  of  defects  including  microcracks,  pores,  and  weak  interfacial  zones. 
Additional  microcracks  nucleate  and  grow  due  to  stresses  induced  by  me¬ 
chanical  or  tectonic  loading,  differential  expansion/contraction  between 
grains  from  thermal  loading,  expansion  of  fluids  in  void  space  during 
freezing/thawing,  or  by  various  chemical  processes.  The  presence  of  mi¬ 
crocracks  in  the  material  affects  the  macroscale  behavior  of  a  material  by 
causing  non-linearity  of  stress-  strain  relations,  reducing  the  stiffness  and 
strength,  inducing  material  anisotropy,  and  causing  irreversible  strains 
due  to  cracks  opening  among  other  effects  (Shao  and  Rudnicki  2000). 
Many  studies  have  shown  the  effects  of  microcracking  on  quasi-brittle 
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material  properties.  This  review  will  focus  on  those  showing  the  effects  of 
microcracking  caused  by  mechanical  and  thermal  loading  because  those 
studies  are  most  relevant  to  this  work. 

Microcrack  properties  cannot  be  used  to  determine  the  stress  and  defor¬ 
mation  history  of  a  material  (Wawersik  and  Brace  1971).  However,  mi¬ 
crocrack  properties  can  be  used  to  predict  the  mechanical  properties  of  a 
material,  and  there  remains  a  need  to  correlate  microscale  phenomena 
with  the  macroscale  mechanical  behavior.  To  meet  that  need,  this  work 
seeks  to  link  the  mechanical  response  of  limestone  to  physical  changes  at 
the  microstructural  level  with  a  view  to  improving  the  predictive  capabili¬ 
ties  of  engineering  material  response  models.  This  includes  providing 
much  needed  mechanical  test  data  for  materials  with  varying  levels  of  mi¬ 
crocracks.  One  of  the  intentions  of  this  work  is  to  establish  experimental 
methods  that  can  be  applied  to  all  quasi-brittle  materials;  as  such,  this  re¬ 
view  will  cover  the  effects  of  microcracking  for  a  variety  of  quasi-static  ma¬ 
terials.  For  more  detailed  information  and  a  review  of  microcracking,  see 
Kranz  (1983). 

2.2  Mechanically  induced  microcracking 

During  the  mechanical  loading  of  quasi-brittle  materials,  microcracks  nu¬ 
cleate,  grow,  and  coalesce  under  deviatoric  stresses  (as  opposed  to  hydro¬ 
static  stresses  that  mainly  affect  pores)  to  cause  failure.  Initially, 
microcracks  increase  in  the  weak  interfacial  zone  between  aggregates  and 
matrix  material  or  along  grain  boundaries.  As  loading  increases,  mi¬ 
crocracks  coalesce  into  continuous  mesoscopic  cracks,  usually  parallel  to 
the  major  principal  axis  in  compression.  Approaching  the  failure  stress, 
the  micro-  and  mesoscale  cracking  leads  to  a  localization  of  cracks  com¬ 
posed  of  a  system  of  bridging  macrocracks.  This  localization  continues  to 
total  failure  by  the  development  of  a  single  macrocrack  (Mertens  2009). 
For  example,  microcracking  during  uniaxial  loading  begins  when  the 
stress  level  is  between  30  and  50  percent  of  the  compressive  strength  in 
compression  (Brace  et  al.  1966)  and  between  50  and  80  percent  of  the  ten¬ 
sile  strength  in  tension  (Ockert  1997  in  Riedel  et  al.  1999).  Above  that 
stress  level,  microcracks  continue  to  develop  and  coalesce  into  a  mac¬ 
roscale  crack  at  failure. 

Studies  have  examined  the  relationship  between  the  mechanical  loading 
process  and  microcracking  during  different  stages  of  mechanical  loading 
using  optical  microscopy  (OM)  and  scanning  electron  microscopy  (SEM). 
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All  of  these  studies  analyzed  microcracks  from  samples  that  had  been  cre¬ 
ated  following  mechanical  loading,  unloading,  and  sample  preparation  for 
the  observational  studies.  The  studies  most  relevant  to  this  work  are  dis¬ 
cussed  next.  The  particular  findings  of  these  works  are  not  the  focus  of  the 
review,  but  rather  the  methods  used  to  obtain  the  microscale  properties  af¬ 
ter  mechanical  loading. 

Hadley  (1976)  examined  microcracks  prior  to  failure  for  Westerly  granite 
samples  loaded  in  triaxial  compression  up  to  80  MPa  confining  pressure. 
The  samples  stressed  from  the  triaxial  tests  were  analyzed  with  SEM  using 
a  direct  measurement  approach. 

The  direct  measurement  approach  consists  of  directly  measuring  crack 
lengths  and  orientations  on  a  plane  section  (Wong  1985)  and  was  used  in 
this  work  to  determine  crack  lengths,  widths,  and  orientations. 

Golshani  et  al.  (2006)  also  used  the  direct  measurement  approach  to  de¬ 
termine  microcrack  lengths  and  densities  of  samples  loaded  in  uniaxial 
and  triaxial  compression.  Confining  pressure  of  20,  40,  60,  and  80  MPa 
were  used  and  microcrack  properties  were  measured  at  various  deviatoric 
stresses  for  each  confining  pressure.  The  microcrack  measurements  were 
used  to  validate  the  material  model  presented  in  the  paper. 

Wawersik  and  Brace  (1971)  quantified  microcracks  from  two  rocks  (West¬ 
erly  granite  and  Frederick  diabase)  that  had  been  loaded  to  different 
stages  in  the  post-failure  region  in  uniaxial  and  triaxial  compression.  A 
mechanical  testing  technique  was  used  to  load  the  samples  at  and  beyond 
the  ultimate  strengths  of  the  materials.  The  confining  pressures  used  went 
up  to  approximately  150  MPa  for  granite  and  15  MPa  for  diabase.  Follow¬ 
ing  the  post-failure  loading,  polished  sections  were  made  for  OM  analysis 
using  a  geometric  probability  approach  for  crack  measurements.  The  geo¬ 
metric  probability  approach  involves  counting  the  number  of  intersections 
between  microcracks  in  a  plane  section  with  an  array  of  parallel  lines 
(Wong  1985).  The  microstructural  analysis  yielded  information  on  crack 
densities,  crack  orientations,  and  their  location  relative  to  grain  bounda¬ 
ries. 

Homand  et  al.  (2000)  examined  microcracks  prior  to  failure  for  Vienne 
granite  samples  loaded  in  triaxial  compression  at  5  MPa  confining  pres¬ 
sure.  OM  analyses  were  done  on  intact  samples  and  samples  loaded  to  60, 
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80,  90,  99,  and  100  percent  of  the  strength  of  the  material.  The  geometric 
probability  approach  was  used  to  measure  crack  densities  and  orientations 
from  these  tests. 

Fredrich  et  al.  (1989)  also  used  the  geometric  probability  approach  to  ex¬ 
amine  pre-failure  samples  of  Carrarra  marble.  The  samples  were  loaded  in 
triaxial  compression  at  various  deviatoric  stresses  and  at  confining  pres¬ 
sures  up  to  450  MPa.  After  samples  had  been  loaded  to  the  desired  stress 
state  and  unloaded,  crack  densities  in  two  orthogonal  directions  were 
measured  using  OM. 

Tapponnier  and  Brace  (1976)  used  the  geometric  probability  approach  to 
examine  samples  triaxially  tested  in  the  pre-  and  post-failure  regions.  For 
the  pre-failure  tests,  samples  were  loaded  to  stresses  at  various  increments 
of  the  peak  strength  of  the  material,  unloaded,  and  observed  using  SEM. 
The  post-failure  samples  were  those  created  by  Wawersik  and  Brace 
(1971). 

Each  of  these  six  studies  analyzed  microcracks  after  loading  to  a  particular 
stress  state,  unloading,  and  sample  preparation  for  the  microscale  anal¬ 
yses.  When  analyzing  samples  created  during  this  process,  it  is  unclear 
whether  the  microscale  properties  taken  after  the  load/unload  cycle  to  a 
particular  stress  represents  the  actual  microscale  properties  at  that  stress 
state.  Before  using  this  data  to  determine  a  correlation  between  the  quan¬ 
tified  microscale  properties  and  its  effects  on  the  mechanical  properties, 
this  method  should  be  validated  to  show  that  the  microscale  properties  de¬ 
termined  from  this  method  represent  the  microstructural  state  of  the 
stress  state  from  which  they  were  taken.  Nonetheless,  mechanical  testing 
can  be  used  to  induce  microcracking  for  the  purpose  of  correlating  micro- 
structural  properties  with  mechanical  response  as  long  as  the  mechanical 
response  is  tested  following  the  microstructural  characterization. 

Zhao  (1998)  used  real-time  SEM  during  loading  to  observe  microcracking 
in  Fangshen  marble.  Microcrack  lengths  and  orientations  were  directly 
measured  during  uniaxial  loading  using  SEM  on  a  20  mm  x  5  mm  x  2  mm 
sample  (4  mm2  area  of  observation).  The  data  from  this  method  is  pre¬ 
ferred  over  that  taken  from  samples  that  have  been  subjected  to  a  load/un¬ 
load  cycle  because  the  microcracking  observed  during  loading  is  clearly 
correlated  to  a  particular  stress  state.  However,  this  method  can  only  be 
used  for  uniaxial  compression  tests  limiting  its  applicability.  This  method 
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is  also  limited  to  small  samples,  and  it  may  be  difficult  to  test  materials 
with  large  aggregates  or  grain  sizes. 

To  link  the  mechanical  response  of  limestone  to  physical  changes  at  the 
microstructural  level,  this  work  induces  various  levels  of  microcracking 
prior  to  mechanical  testing.  This  method  avoids  the  uncertainty  associated 
with  the  results  obtained  following  a  load/unload  cycle.  This  method  can 
also  be  used  on  relatively  large  samples  and  to  test  samples  under  multi- 
axial  stress  states,  which  cannot  be  done  using  real-time  observation  tech¬ 
niques.  Though  mechanical  testing  is  a  viable  option  to  induce  microcrack¬ 
ing,  thermal  treating  is  used  to  induce  various  levels  of  microcracking  in 
this  work  for  reasons  mentioned  later. 

2.3  Thermally  induced  microcracking 

When  quasi-brittle  materials  are  subject  to  thermal  loading  (i.e.,  heating 
to  a  maximum  temperature  followed  by  cooling  back  to  room  tempera¬ 
ture),  microcracking  is  induced  by  slight  differences  in  the  coefficient  of 
thermal  expansion  between  grains,  which  cause  differential  expansion  or 
contraction  under  sufficient  thermal  load.  In  rocks,  new  microcracks  are 
formed  after  the  previous  maximum  temperature  has  been  exceeded 
(Johnson  et  al.  1978;  Yong  and  Wang  1980).  Slow,  uniform  heating  above 
the  previous  maximum  temperature  cause  grain  boundary  separation  and 
intragranular  cracking  (Sprunt  and  Brace  1974;  Simmons  and  Richter 
1976;  Friedman  and  Johnson  1978;  Bauer  and  Johnson  1979;  Kranz  1983) 
while  large  heating  rates  result  in  thermal  gradients  that  can  cause  addi¬ 
tional  microcracking  (Johnson  and  Gangi  1980). 

To  determine  the  effects  of  thermal  loading  on  microstructural  properties, 
studies  have  been  performed  to  investigate  microcracking  using  OM  and 
SEM  following  thermal  treating  to  various  maximum  temperatures.  Some 
studies  (Fredrich  and  Wong  1986;  Lin  2002)  quantified  microcracks  after 
various  thermal  treatments  but  did  not  perform  mechanical  tests  to  deter¬ 
mine  a  correlation  between  the  microcracking  and  changes  in  mechanical 
properties.  Other  studies  (Friedman  and  Johnson  1978;  Homand-Etienne 
and  Troalen  1984;  Yavuz  et  al.  2010)  provided  only  qualitative  changes  in 
the  materials  following  heat  treatments  and  also  did  not  perform  mechani¬ 
cal  tests.  Because  this  work  is  investigating  the  correlation  between  mi¬ 
crocracking  and  its  effects  on  the  mechanical  response,  only  those  studies 
that  investigated  microcracks  with  corresponding  mechanical  tests  will  be 
reviewed. 
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Keshavarz  et  al.  (2010)  performed  uniaxial  compression  tests  on  gabbro 
samples  heated  up  to  1000  °C.  The  uniaxial  compression  data  showed  a 
steady  decrease  in  strengths  up  to  6oo°C  and  a  sharp  decrease  thereafter. 
The  qualitative  OM  analysis  showed  that  microcracks  were  the  main 
source  of  strength  reduction  up  to  6oo°C,  and  phase  transformations  (in 
this  case,  oxidation  of  certain  minerals)  and  bursting  of  fluid  inclusions 
(there  was  no  indication  that  materials  were  dried  before  heating)  caused 
an  increased  strength  reduction  rate  thereafter.  However,  microcracking 
was  not  quantified  in  this  study  so  no  correlation  can  be  made  between  the 
microstructural  changes  and  its  effects  on  the  mechanical  properties. 

Ferrero  and  Marini  (2001)  performed  uniaxial  compression  tests  on  two 
marbles  (Ormea  black  marble  and  Perlato  Sicilia  marble)  heated  up  to 
6oo°C.  The  mechanical  test  data  showed  decreasing  elastic  moduli,  tensile 
strengths,  and  sonic  velocities  with  increasing  maximum  temperatures  in 
the  heat  treatments.  The  OM  analysis  on  each  material  also  showed  an  in¬ 
crease  in  microcrack  density  for  each  material  with  increasing  tempera¬ 
tures. 

Homand-Etienne  and  Houpert  (1989)  performed  uniaxial  compression, 
hydrostatic  compression,  and  uniaxial  tension  tests  on  two  granites 
(Senones  and  Remiremont  granite)  heated  up  to  6oo°C.  The  study  showed 
that  crack  density  increases  and  elastic  moduli  and  strengths  decrease  as 
function  of  the  maximum  temperatures  obtained  in  the  heat  treatments. 
The  study  also  found  the  elastic  moduli  to  be  more  sensitive  to  mi- 
crocracks  than  the  compressive  strengths. 

Lastly,  Bauer  and  Johnson  (1979)  performed  uniaxial  and  triaxial  com¬ 
pression  (confining  pressures  up  to  60  MPa)  and  uniaxial  tension  tests  on 
two  granites  (Westerly  and  Charcoal  granite).  The  study  found  that  below 
500°C  formation  of  new  microcracks  was  the  primary  microstructural 
change;  above  500°C,  cracks  widened  and  the  material  went  through  a 
phase  transformation.  Regarding  the  mechanical  testing,  the  tensile 
strengths  of  the  granites  did  not  decrease  until  maximum  heat  treatment 
temperatures  exceeded  200°C.  The  compression  tests  showed  that  confin¬ 
ing  pressure  reduced  the  effects  of  microcracks  on  the  compressive 
strength. 

Each  of  these  studies  showed  strength  and  elastic  moduli  decreases  with 
increasing  temperatures.  Bauer  and  Johnson  (1979)  also  observed  that  as 


ERDC/GSL  TR-17-15 


10 


confining  pressures  increased,  the  effects  of  microcracking  on  strength  de¬ 
creased.  Bauer  and  Johnson  (1979),  Homand-Etienne  and  Houpert 
(1989),  and  Ferrero  and  Marini  (2001)  performed  quantitative  microcrack 
analyses  along  with  their  mechanical  tests,  so  these  studies  could  be  used 
to  correlate  microcracking  with  macroscale  material  behavior. 

Similar  to  these  studies,  this  work  characterizes  microcracking  in  Salem 
limestone  following  thermal  treating  to  various  maximum  temperatures. 
However,  the  focus  of  this  work  is  to  provide  data  on  a  material  with  vary¬ 
ing  levels  of  microcracks  (and  methods  to  test  other  quasi-brittle  materi¬ 
als)  with  a  view  to  improving  predictive  capabilities  of  engineering 
material  response  models.  With  this  focus  in  mind,  the  thermal  treatment 
was  limited  to  maximum  temperatures  below  which  microcracking  is  the 
only  microscale  phenomenon  occurring.  This  ensures  that  the  quantified 
microcrack  data  and  corresponding  mechanical  tests  could  also  be  corre¬ 
lated  with  mechanically  induced  microcracking.  This  work  also  provides 
additional  mechanical  tests  than  previous  studies;  the  mechanical  tests 
done  in  this  work  includes  uniaxial  compression  (quasi-static  and  high 
strain  rate),  triaxial  compression  (confining  pressures  up  to  400  MPa),  hy¬ 
drostatic  compression,  and  uniaxial  strain/constant  volume  strain  path 
tests  (explained  in  later  chapters).  This  data  set  can  be  used  to  provide  in¬ 
sight  on  quasi-brittle  material  behavior  at  multiple  stress  states  and  strain 
rates,  which  is  vitally  important  when  modeling  these  materials. 

2.4  Quasi-brittle  material  modeling 

The  ability  to  model  the  behavior  of  structures  consisting  of  quasi-brittle 
materials,  such  as  concrete  and  rock,  is  important  in  many  civil  engineer¬ 
ing  applications  due  to  their  wide  use  in  civil  and  military  structures. 

These  models  are  necessary  for  the  entire  life  cycle  of  the  structure  begin¬ 
ning  at  the  design  stage  and  extending  throughout  the  life  of  structure,  in¬ 
cluding  renovations  and  repairs  following  damaging  events. 

Quasi-brittle  materials  exhibit  complex  behavior  that  is  often  non-linear 
and  changes  with  varying  stress  states,  strain  rates,  and  damage.  Each  of 
these  aspects  must  be  taken  into  account  in  order  to  adequately  design 
structures  made  from  quasi-brittle  materials. 

Numerical  models  generally  used  to  predict  the  behavior  of  quasi-brittle 
materials  are  categorized  as  either  micromechanical  models  or  phenome¬ 
nological  models. 
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Micromechanical  models  attempt  to  model  real  physical  mechanisms  in¬ 
volved  in  inelastic  deformation,  such  as  microcrack  nucleation  and  propa¬ 
gation,  but  are  not  suited  for  structural  scale  problems  because 
representative  volumes  of  the  material  are  explicitly  modeled  (Hoxha  and 
Homand  2000).  Phenomenological  models  are  developed  in  the  irreversi¬ 
ble  thermodynamic  framework  with  relatively  simple  macroscopic  consti¬ 
tutive  equations  (Shao  et  al.  1999;  Hoxha  and  Homand  2000;  Shao  and 
Rudnicki  2000)  and  do  take  into  account  the  origin  of  the  damage.  These 
models  comprise  many  different  types  of  models  including  classical  plas¬ 
ticity  models,  hypoelastic  models,  nonlinear  elastic  models,  viscoplastic 
models,  and  continuum  damage  mechanics  models  among  others  (Bazant 
and  Prat  1988). 

Because  this  work  is  intended  for  use  on  the  structural  scale,  this  review 
will  look  at  two  types  of  phenomenological  models,  continuum  damage 
mechanics  models  and  plasticity  models.  The  intent  of  this  section  is  not  to 
provide  an  exhaustive  review  of  all  phenomenological  models  but  to  pro¬ 
vide  an  idea  of  how  these  two  types  of  models  handle  damage  and  how 
they  can  be  improved  by  the  experimental  data  presented  in  this  work. 

To  show  how  quantified  microcrack  data  is  used  in  brittle  material  models, 
the  continuum  damage  mechanics  models  of  Hoxha  and  Homand  (2000), 
Shao  and  Rudnicki  (2000),  and  Golshani  et  al.  (2006)  will  be  discussed 
first.  Due  to  the  inherent  complexity  of  continuum  damage  mechanics 
models  and  the  experimental  focus  of  this  work,  a  general  description  of 
these  models  is  given  without  their  theoretical  basis.  Then,  the  commonly 
used  Holmquist-Johnson-Cook  concrete  model  (Holmquist  et  al.  1993) 
will  be  discussed,  which  will  give  insight  into  damage  formulations  in  a 
plasticity  model. 

2.4.1  Continuum  damage  mechanics  models 

In  continuum  damage  mechanics  models,  microscopic  events  during  load¬ 
ing  such  as  microcrack  initiation,  growth,  and  kinetics  are  represented  in 
the  macroscopic  response  by  scalar  or  tensor  variables.  Because  the  micro¬ 
scopic  events  are  taken  into  account,  these  models  can  offer  information 
on  microcrack  related  quantities,  e.g.,  microcrack  density,  length,  and  ori¬ 
entation,  during  loading  (Golshani  et  al.  2006).  However,  it  is  difficult  to 
quantify  the  microcracking  and  apply  it  to  applications  using  these  models 
(Zhou  2010).  Of  this  type  of  model,  three  of  the  most  relevant  models  to 
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this  work  will  be  discussed  further,  Hoxha  and  Homand  (2000),  Shao  and 
Rudnicki  (2000),  and  Golshani  et  al.  (2006). 

Hoxha  and  Homand  (2000)  proposed  a  micromechanics-based  damage 
model  that  uses  a  stereological  parameter  to  capture  the  changes  in  crack 
geometry.  The  stereological  parameter  is  the  number  per  unit  length  of 
crack  intersections  with  a  test  line  defined  for  a  given  direction.  They  used 
the  published  results  of  previously  mentioned  microscopic  analyses  (Had¬ 
ley  1976;  Wong  1985;  Zhao  1998;  Homand  et  al.  2000)  to  formulate  their 
microcrack  evolution  equations.  The  evolution  of  the  crack  geometry  is 
used  to  predict  the  effective  properties  of  the  material  during  compressive 
loading.  This  model  can  provide  good  agreement  with  laboratory  tests  un¬ 
til  the  last  stage  of  testing  when  crack  coalescence  begins.  At  that  point, 
the  simulations  deviate  from  the  test  data.  This  model  was  only  used  for 
compressive  stress  states. 

Shao  and  Rudnicki  (2000)  have  also  proposed  a  micromechanics-based 
damage  model,  which  uses  an  internal  variable  to  represent  the  density 
and  orientation  of  microcracks.  They  performed  their  own  tests  to  induce 
damage  for  their  proposed  microcrack  evolution  that  consisted  of  confined 
compression  tests  at  different  confining  pressures.  This  macroscale  data 
and  various  assumptions  about  the  microcrack  growth  are  used  to  corre¬ 
late  the  macroscale  data  with  microscale  events.  No  efforts  were  made  to 
verify  whether  the  assumed  microscale  events  were  actually  occurring.  As 
with  Hoxha  and  Homand  (2000),  the  model  works  well  until  the  samples 
are  near  failure.  At  this  point,  the  model  is  not  in  agreement  with  experi¬ 
mental  data. 

The  proposed  model  of  Golshani  et  al.  (2006)  predicts  changes  in  mi¬ 
crocrack  lengths  and  density  and  compares  those  microcrack  changes  with 
experimental  data. 

They  performed  their  own  tests  to  induce  and  quantify  damage,  which 
consisted  of  confined  compression  tests  at  different  confining  pressures. 
During  mechanical  testing,  the  deviatoric  stresses  were  varied  in  order  to 
monitor  crack  property  changes  at  the  various  deviatoric  stresses.  As  with 
the  other  two  models,  the  final  stages  of  testing  do  not  compare  well  with 
experimental  data,  but  the  model  was  able  to  predict  the  overall  trend  of 
changes  in  microcrack  length  and  density  under  triaxial  loading. 
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The  models  by  Hoxha  and  Homand  (2000)  and  Shao  and  Rudnicki  (2000) 
account  for  measurable  microcrack  properties  and  their  evolution,  but 
they  do  not  provide  evidence  that  the  models  predict  the  actual  microcrack 
evolution  during  loading.  Only  the  effects  on  the  macroscale  response  are 
given  in  the  form  of  stress-strain  relationships.  The  model  of  Golshani  et 
al.  (2006)  showed  how  microcracks  evolved  with  stress  along  with  the 
macroscale  stress-strain  relationships.  Golshani  stated,  ‘any  microme¬ 
chanical  model  is  not  fully  justified  even  if  the  predictions,  on  a  macro¬ 
scale,  agree  well  with  experimental  results.’  Though  the  initial  damage 
state  may  be  correctly  implemented,  the  damage  evolution  should  be  com¬ 
pared  with  experimental  results  to  ensure  the  model  is  capturing  the  dam¬ 
age  evolution  correctly. 

2.4.2  Plasticity  models 

The  complexity  of  continuum  damage  mechanics  has  mostly  limited  their 
use  to  research  endeavors.  In  constrast,  the  relative  simplicity  of  plasticity 
models  has  resulted  in  more  widespread  use.  The  simplicity  of  the  plastic¬ 
ity  models  comes  from  their  implementation  of  relatively  simple  macro¬ 
scopic  constitutive  equations  that  do  not  take  into  account  the  origin  of 
damage.  To  show  how  the  experimental  data  obtained  in  this  work  can  be 
used  to  improve  material  models  in  general,  the  Holmquist-Johnson-Cook 
(HJC)  model  (1993),  a  popular  material  model  for  impact  problems,  has 
been  chosen  for  modification  based  on  the  material  response  observed  in 
this  work. 

The  HJC  model  is  an  elastic- viscoplastic  model  with  isotropic  damage  for 
concrete  subjected  to  large  strains  and  high  strain  rates  and  pressures. 

This  model  is  included  because  it  has  been  widely  used,  modified  (Riedel 
et  al.  1999;  Gebbeken  and  Ruppert  2000;  Polanco-Loria  et  al.  2008),  and 
takes  into  account  the  most  important  issues  of  brittle  material  behavior, 
such  as  pressure  dependency,  rate  dependency,  pressure-volume  changes, 
and  damage. 

In  the  HJC  model,  the  material  is  linear  elastic  until  the  failure  surface,  or 
shear  limit  surface,  is  reached  and  damage  is  accumulated  by  inelastic  de¬ 
formations.  Damage  then  affects  the  material  by  reducing  the  cohesive 
strength  value,  A,  in  the  failure  surface  given  in  Eqs.  2.1  and  2.2. 


ERDC/GSL  TR-17-15 


14 


a 


[A(l  -  D)  +  BP*W]  [1  +  C  Inf*]  for  a*  <  SMAX 


(2.1) 


a*  =  SMAX  for  a*  >  SMAX 


(2.2) 


In  Eqs.  2.1  and  2.2,  o  *  and  P*  are  the  deviatoric  stress  and  pressure  nor¬ 
malized  with  respect  to  the  unconfined  compressive  strength  of  the  mate¬ 
rial, /''c  ,  respectively,  and  e*  is  the  strain  rate  normalized  to  a  reference 
strain  rate  (usually  1  s_1).  A,  B,  N,  and  C  are  material  constants,  and  SMAX 
is  the  maximum  strength.  D  represents  damage  and  is  given  by 


q  _  y  (Afp+AMp) 


(2.3) 


where  AfP  and  ApP  are  the  equivalent  plastic  strain  increment  and  equiva¬ 
lent  volumetric  strain  increment,  respectively,  and  efp  and  /ifp  those  strains 

at  failure.  The  failure  surface  and  the  effect  of  damage  on  it  are  shown  in 
Figure  2-1. 


Figure  2-1.  The  shear  limit  surface  of  the  HJC  model. 


* 


P* 


Though  this  model  takes  into  account  some  important  aspects  of  concrete 
behavior,  some  behavior  is  not  accounted  for  such  as  the  non-linearity  be¬ 
fore  peak  load,  stiffness  degradation  and  recovery,  induced  anisotropy, 
and  shear  strength  differences  between  tension  and  compression  (Po- 
lanco-Loria  et  al.  2008).  There  is  no  evidence  to  support  that  the  damage 
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effects  implemented  in  the  HJC  model,  given  in  Eqs.  l  and  2  and  shown  in 
Figure  1,  reflect  what  actually  occurs  when  the  material  is  damaged. 

To  improve  the  implementation  of  damage  in  the  HJC  model,  basic  princi¬ 
ples  determined  from  the  experimental  data  performed  in  this  work  were 
used  to  make  modifications  to  the  model.  The  goal  of  these  modifications 
was  to  improve  the  prediction  capabilities  for  damaged  materials  under 
varying  stress  states  and  strain  rates.  For  comparison,  results  from  the 
original  HJC  model  and  the  modified  version  are  compared  to  the  experi¬ 
mental  data. 

2.5  Summary 

The  first  portion  of  this  review  covered  microcracking,  the  primary  dam¬ 
age  mechanism  in  quasi-brittle  materials.  The  effects  of  both  mechanically 
and  thermally  induced  microcracking  were  discussed,  which  included 
studies  that  correlated  quantified  microcracks  with  changes  in  macroscale 
mechanical  behavior.  The  studies  using  mechanical  loading  to  induce  mi¬ 
crocracking  performed  their  microstructural  analysis  after  a  loading/un- 
loading  cycle,  which  leads  to  some  uncertainty  about  how  the 
microstructural  characterization  is  related  to  the  stress  state  from  which  it 
was  taken. 

Regarding  the  studies  on  thermally  induced  microcracking,  the  effects  of 
microcracking  on  a  limited  number  of  stress  states  were  investigated. 

To  address  these  issues,  this  work  induced  various  levels  of  microcracking 
prior  to  mechanical  testing  that  included  uniaxial  compression  (quasi¬ 
static  and  high  strain  rate),  triaxial  compression  (confining  pressures  up  to 
400  MPa),  hydrostatic  compression,  and  uniaxial  strain/ constant  volume 
strain  path  tests.  This  data  set  can  be  used  to  provide  insight  on  micro- 
and  macroscale  quasi-brittle  material  behavior  at  multiple  stress  states 
and  strain  rates,  which  is  necessary  when  modeling  these  materials. 

Four  phenomenological  material  models  have  been  discussed  including 
three  continuum  damage  mechanics  models  and  a  plasticity  model.  The 
recent  development  of  micromechanical-based  damage  models  have  al¬ 
lowed  for  internal  variables  based  on  physical  mechanisms  such  as  mi¬ 
crocrack  densities  and  orientations.  However,  there  is  still  much  work 
before  these  can  be  used  in  a  predictive  capacity.  They  work  well  for  pre¬ 
dicting  overall  effective  stiffness  but  cannot  predict  failure  because  it  is 
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usually  due  to  the  growth  of  cracks  of  one  orientation  (Hoxha  and  Ho- 
mand  2000)  or  frictional-dilatant  slip  in  a  particular  direction  (Bazant  and 
Zi  2003).  The  HJC  plasticity  model  was  also  discussed  because  of  its  wide¬ 
spread  use  and  simplicity  of  its  constitutive  equations. 

Though  the  HJC  model  captures  important  aspects  of  concrete  behavior, 
there  are  still  other  important  aspects  that  should  be  considered.  Strength 
and  stiffness  degradation  due  to  damage  are  the  aspects  addressed  in  this 
work. 

There  is  a  very  limited  amount  of  experimental  data  with  which  to  cali¬ 
brate  and  validate  any  of  the  material  models  mentioned.  The  micro-  and 
macroscale  data  provided  by  this  work  will  aid  in  the  development  of  con¬ 
tinuum  damage  mechanics  models  by  providing  additional  data  sets  for 
use  in  model  validation  and  calibration.  The  data  in  this  work  should  also 
aid  in  the  development  of  plasticity  models  by  providing  basic  principles 
with  which  to  formulate  the  macroscale  response  of  a  damaged  material. 
To  show  how  this  work  can  be  used  to  improve  phenomenological  plastic¬ 
ity  models,  the  effects  of  damage  characterized  in  this  work  are  used  to 
modify  the  damage  implementation  in  the  HJC  concrete  plasticity  model. 
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3  Inducing  and  Quantifying  Damage 

3.1  Overview 

This  chapter  is  divided  into  three  sections:  Indiana  limestone  description; 
the  experimental  method  used  to  induce  damage;  and  the  experimental 
methods  used  to  quantify  the  induced  damage.  After  the  material  descrip¬ 
tion,  the  thermal  treating  method  used  to  induce  damage  in  the  samples 
will  be  discussed.  Then,  the  methods  used  to  quantify  that  damage  will  be 
given  along  with  the  results  of  the  quantification. 

3.2  Materials 

The  material  used  in  this  study  is  commonly  referred  to  as  Indiana,  Salem, 
or  Bedford  limestone.  The  material  is  of  Mississippian  age  and  is  a  com¬ 
monly  used  building  stone  material.  It  was  quarried  and  cut  by  the  Elliot 
Stone  Company  of  Bedford,  IN,  into  i8-in.  by  i8-in.  by  12-in.  blocks  (see 
Figure  3-1).  Indiana  limestone  is  composed  primarily  of  calcium  carbonate 
(CaC03),  or  calcite,  that  forms  the  grains  and  the  cementing  material  that 
binds  the  grains.  The  limestone  is  uniform  in  composition,  texture,  and 
structure  and  is  practically  non-crystalline  in  character  (ILI  Handbook 
2007). 


Figure  3-1.  Salem  limestone  blocks  from  Elliot  Stone  Company. 


The  minimum  compressive  strength  of  Indiana  limestone  (ILI  Handbook 
2007)  is  approximately  27.6  MPa  (4000  psi).  Frew  et  al.  (2010)  found  the 
unconfined  compressive  strength  to  be  approximately  65  MPa  (9400  psi). 
Porosity  varies  between  12  and  21  percent,  depending  on  quarry  location 
(Churcher  et  al.  1991;  Fossum  et  al.  1995).  The  ILI  Handbook  (2007)  also 
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gives  a  range  for  bulk  specific  gravity  between  2.1  and  2.75,  and  modulus 
of  elasticity  between  22,750  and  37250  MPa  (3,300,000-5,400,000  psi). 
The  material  used  in  this  work  had  an  unconfined  compressive  strength  of 
72  MPa  (10400  psi). 

This  particular  material  was  chosen  for  this  work  for  three  reasons:  (1) 
There  have  been  several  other  studies  done  on  Salem  limestone  that  pro¬ 
vide  additional  test  data  (Cummins  1991;  Frew  2001;  Frew  et  al.  2010);  (2) 
Its  mechanical  behavior  should  be  analogous  to  other  quasi-brittle  materi¬ 
als,  e.g.,  concrete,  ceramics,  and  other  rocks;  and 

(3)  This  limestone  is  relatively  homogeneous  compared  to  other  quasi- 
brittle  materials,  which  was  desired  when  developing  the  test  method  used 
in  this  work. 

3.3  Inducing  damage 

In  order  to  study  the  effects  of  microcracking  on  the  mechanical  properties 
of  limestone,  a  method  of  inducing  consistent  levels  of  damage  in  the  lime¬ 
stone  samples  is  needed.  Mechanical  loading,  thermal  treating,  freez¬ 
ing/thawing,  and  chemical  treating  have  been  used  to  investigate  the 
processes  involved  in  mechanical  or  tectonic  loading,  thermal  loading, 
weathering,  and  corrosion,  respectively.  However,  thermal  treating  was 
chosen  for  the  reasons  discussed  next. 

Thermal  treating  is  the  process  of  heating  a  material  to  a  maximum  tem¬ 
perature  and  cooling  the  sample  back  to  room  temperature.  This  process 
induces  microcracks  because  slight  differences  in  the  coefficient  of  thermal 
expansion  between  grains  cause  local  stresses.  If  the  maximum  tempera¬ 
ture  is  kept  below  that  which  causes  changes  in  the  material  composition, 
the  microstructural  changes  caused  by  the  heat  treatments  should  be  com¬ 
parable  with  that  caused  by  mechanical  loading  and  freezing/thawing. 

Chemical  treating  of  the  samples  was  not  chosen  to  induce  microcracking 
because  it  may  cause  changes  that  are  not  related  to  the  damage  mecha¬ 
nisms  found  in  mechanical  loading,  thermal  treating,  and  freezing/thaw- 
ing. 

Thermal  treating  was  chosen  to  induce  microcracking  rather  than  mechan¬ 
ical  loading  and  freezing/thawing  because  it  is  less  complicated  and  less 
expensive  to  perform.  When  using  thermal  treatment,  heating/cooling  rate 
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and  maximum  temperature  need  to  be  considered,  and  a  low  heating/cool¬ 
ing  rate  is  usually  chosen  to  reduce  thermal  shock  (Homand-Etienne  and 
Houpert  1989;  Homand-Etienne  and  Troalen  1984).  If  freezing/thawing  is 
used  to  induce  damage,  the  degree  of  saturation,  freezing/thawing  rate, 
and  minimum  temperature  are  the  major  variables  concerned.  Constant 
saturation  throughout  the  sample  also  needs  to  be  carefully  handled  or  lo¬ 
calized  damage  could  occur.  Regarding  the  costs  associated  with  thermal 
treating,  an  oven  with  controllable  heating  cooling/rate  is  all  that  is 
needed,  and  the  samples  can  be  treated  simultaneously  in  about  24  hours. 
For  mechanical  testing  (e.g.,  triaxial  testing  at  400  MPa  confining  pres¬ 
sures),  two  samples  can  be  tested  a  day  for  approximately  $1500  in  labor 
alone. 

Lastly,  thermal  treating  can  be  used  to  induce  various  levels  of  microcrack¬ 
ing  by  simply  varying  the  maximum  temperature  reached  during  the  ther¬ 
mal  treatment.  To  get  consistent  levels  of  damage  while  running 
mechanical  tests,  the  sample  strengths  must  be  fairly  reproducible  and  the 
samples  must  be  closely  monitored  when  nearing  peak  stress.  Most  me¬ 
chanical  tests  also  introduce  anisotropic  microcracking,  whereas  heat 
treating  induces  relatively  isotropic  microcracking  (David  et  al.  1999). 

3.3.1  Thermal  treating 

To  induce  microcracks  in  the  limestone,  the  samples  were  heated  in  a  muf¬ 
fle  furnace  at  a  low  heating  rate  (i°C/min  to  300°C,  2°C/min  thereafter). 
The  maximum  temperature  was  then  held  for  five  hours  and  cooled  down 
at  approximately  1  °C/min.  This  heat  treatment  was  chosen  because  it  is 
similar  to  heat  treatments  used  in  previous  studies  (Homand-Etienne  and 
Houpert  1989;  Homand-Etienne  and  Troalen  1984).  All  samples  were  kept 
in  a  60°  C  oven  before  and  after  all  mechanical  tests  and  heat  treatments  to 
reduce  their  moisture  content.  Figure  3-2  shows  the  temperature  vs.  time 
for  the  250  °C  and  450  °C  heat  treatments. 
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Figure  3-2.  Temperature  vs.  time  for  the  heat  treatments  used  for  the  mechanical 

test  samples. 


Sample  lengths,  diameters,  and  weights  were  measured  before  and  after 
each  heat  treatment  to  determine  volume  and  density  changes  as  a  func¬ 
tion  of  the  maximum  temperatures  used  in  the  heat  treatments  (see  Figure 
3-3).  These  measurements  were  taken  as  a  rough  estimate  of  the  induced 
microcrack  volume  and  to  give  an  initial  indication  of  other  changes  that 
maybe  occurring  besides  microcracking.  For  the  250  °C  and  450  °C  heat 
treatments,  the  average  percent  volume  changes  were  0.06  and  1.24  per¬ 
cent,  and  the  average  percent  density  changes  were  -0.17  and  -1.48  per¬ 
cent,  respectively.  The  800  °C  heat  treatment  resulted  in  a  drastic  increase 
in  volume  and  decrease  in  density,  which  indicated  that  changes  other 
than  microcracking  were  occurring  around  that  maximum  temperature. 


Figure  3-3.  Percent  volume  change  (left)  and  percent  density  change  vs.  maximum 

heat  treatment  temperature. 
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3.3.2  Wave  velocities 

P-wave  (or  compressional  wave)  velocities  and  S-wave  (or  shear  wave)  ve¬ 
locities  were  measured  using  the  ultrasonic  pulse  transmission  technique 
before  and  after  heat  treatments.  This  technique  is  a  non-destructive 
means  of  evaluating  damage  in  a  material  because  the  waves  passed 
through  the  material  are  influenced  by  the  elastic  properties  and  disconti¬ 
nuities  such  as  pores  and  microcracks.  The  velocities  were  taken  to  give  an 
initial  indication  of  the  effects  of  microcracks  induced  by  the  heat  treat¬ 
ments  and  to  determine  if  heat  treatments  reliably  induce  a  consistent 
amount  of  microcracks. 

For  the  sonic  velocity  tests,  limestone  samples  were  created  using  a  con¬ 
ventional  rock  coring  barrel  to  core  the  original  blocks.  Then,  the  samples 
were  cut  to  their  approximate  final  length  with  a  slabbing  saw  and  later 
ground  flat.  The  final  samples  were  right  cylinders  with  a  length  and  a  di¬ 
ameter  of  50.8  mm  (2.0  in.).  The  experimental  setup  included  a  pair  of  1- 
MHz  piezoelectric  transducers  to  transmit  and  receive  the  P-waves,  a  pair 
of  2.25  MHz  transducers  for  the  S-waves,  and  a  100  MHz  oscilloscope  to 
measure  the  wave  travel  times.  Velocities  were  calculated  by  dividing  the 
length  of  the  samples,  measured  from  digital  micrometers,  by  the  wave 
travel  time.  The  P-  and  S-wave  velocities  were  measured  along  three  or¬ 
thogonal  directions  in  the  cylindrical  samples:  one  in  the  axial  direction 
and  two  in  orthogonal  radial  directions.  This  was  done  to  check  for  uni¬ 
formity  throughout  the  samples.  All  of  the  velocity  measurements  were 
taken  under  atmospheric  conditions  in  accordance  with  ASTM  C  597 
(ASTM  2005c). 

Figure  3-4.  Sonic  velocities  for  increasing  maximum  temperatures. 
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Wave  velocity  measurements,  a  total  of  142,  were  taken  on  the  untreated, 
200  °C,  250  °C,  300  °C,  400  °C,  450  °C,  500  °C,  and  800  °C  heat-treated 
limestone  samples.  The  velocity  data  in  the  axial  direction  are  shown  in 
Figure  3.3.  The  data  show  a  clear  correlation  between  the  heat-treatment 
temperatures  and  the  reduced  wave  velocities  caused  by  induced  mi¬ 
crocracking.  P-wave  velocities  of  untreated  samples  ranged  from  approxi¬ 
mately  4200  m/s  to  4700  m/s,  and  S-wave  velocities  of  untreated  samples 
ranged  from  approximately  2500  m/s  to  2700  m/s.  P-wave  velocities  were 
reduced  by  approximately  17%,  47%,  and  69%,  at  250  °C,  450 °C,  and 
800  °C,  respectively;  and  S-  wave  velocities  were  reduced  by  approximately 
13%,  38%,  and  62%,  at  250 °C,  450  °C,  and  800  °C,  respectively. 

Wave  velocity  measurements  can  also  be  used  to  approximately  calculate 
elastic  moduli.  The  bulk  modulus  and  shear  modulus  were  determined 
from  wave  velocity  measurements  on  undamaged  and  damaged  samples 
using  Eqs.  3.1  and  3.2. 


(3.1) 


G  =  PVS2 


(3.2) 


VP  and  Vs  are  the  P-wave  and  S-wave  velocities  of  the  material,  respec- 
tively.Table  3-1  shows  the  elastic  moduli  for  the  untreated  and  treated 
samples  calculated  from  the  wave  velocities. 

Table  3-1.  Elastic  moduli  for  samples  with  different  heat  treatments. 


Maximum 
Temperature,  °C 

Average  Bulk 
Modulus,  MPa 

Average  Shear 
Modulus,  MPa 

60 

10500 

16000 

200 

7500 

13400 

250 

6500 

12000 

300 

5800 

10700 

400 

3600 

8500 

450 

2100 

6000 

500 

1300 

5000 

800 

600 

1700 

Because  these  values  are  calculated  directly  from  the  wave  velocity  values, 
they  follow  the  same  trend  as  the  wave  velocity,  which  is  a  steady  decrease 
in  elastic  moduli  as  maximum  thermal  treatment  temperatures  increase. 
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These  results  verified  the  initial  assumptions  made  when  choosing  heat 
treatments  for  this  work:  The  heat  treatments  induced  microcracking  as  a 
function  of  the  maximum  temperature,  and  consistent  levels  of  induced 
microcracking  can  be  obtained. 

3.3.3  X-ray  diffraction,  thermogravimetric  analysis,  and  differential 
scanning  calorimetry 

To  exclude  the  possibility  that  the  effects  of  the  heat  treatments  on  the  me¬ 
chanical  response  were  due  to  changes  in  material  composition  (e.g., 
phase  changes),  specimens  were  examined  using  x-ray  diffraction  (XRD), 
thermogravimetric  analysis  (TGA),  and  differential  scanning  calorimetry 
(DSC).  These  tests  were  performed  to  ensure  no  changes  were  occurring 
for  the  temperature  range  of  interest,  ambient  temperature  to  500  °C. 

XRD  is  a  nondestructive  technique  of  determining  the  chemical  composi¬ 
tion  and  crystallographic  structure  of  a  material.  This  method  was  used  to 
detect  material  composition  changes  in  ground  powder  samples  of  the  un¬ 
treated,  200°C,  300°C,  400°C,  500°C,  and  8oo°C  heat-treated  limestone 
samples.  The  XRD  analysis  was  performed  using  a  Panalytical  X’Pert  Pro 
Materials  Research  Diffractometer. 

Figure  3-5  shows  the  XRD  patterns  from  the  limestone  samples  compared 
to  reference  patterns  for  calcite  (CaC03)  and  portlandite  (Ca(0H)2).  The 
8oo°C  heat  treatment  resulted  in  calcite  that  was  partially  transformed  to 
portlandite.  During  the  8oo°C  heat  treatment,  the  calcite  went  through  a 
thermal  decomposition  to  form  calcium  oxide  according  to  the  reaction  in 
Eq.  3.3. 


CaC03  ->  CaO  +  C02  (3.3) 

Upon  removal  of  the  samples  from  the  oven  and  before  the  XRD  test,  the 
calcium  oxide  reacted  with  water  in  the  atmosphere  to  form  portlandite 
according  the  reaction  in  Eq.  3.4. 

CaO  +  H20  ->  Ca(0H)2  (3.4) 

Samples  undergoing  heat  treatments  of  500°C  or  less  show  no  phase 
changes  or  decomposition  of  the  calcite,  which  provides  evidence  that  mi¬ 
crocracking  is  the  main  microstructural  change  for  that  temperature 
range. 
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Figure  3-5.  XRD  results. 
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TGA  is  another  technique  to  detect  material  composition  changes.  This 
method  allows  for  the  determination  of  phases  present  in  the  material 
based  on  the  mass  loss  over  a  range  of  temperatures.  For  this  test,  an  un¬ 
treated  sample  was  analyzed  in  a  nitrogen  atmosphere  (to  prevent  oxida¬ 
tion  and  hydration)  during  a  500°C  heat  treatment  using  the  same  heating 
rates  and  hold  times  for  all  tests  (i°C/min  to  300°C,  and  2°C/min  thereaf¬ 
ter  with  a  5hr  hold  time  at  max  temp).  This  test  was  performed  with  a  Ne- 
tzsch  Jupiter  Thermal  Analyzer.  Figure  3-6  shows  the  results  of  this  test. 

Figure  3-6.  TGA  results. 
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Negligible  changes  in  mass  were  observed,  indicating  no  phase  changes  or 
decomposition  at  temperatures  up  to  500°C.  Significant  changes  in  mate¬ 
rial  composition,  as  detected  by  TGA,  occur  at  temperatures  between 
700°C  and  iooo°C  and  can  cause  mass  losses  of  approximately  45  percent. 

DSC  was  the  last  test  performed  to  detect  phase  changes  in  the  material. 
Material  composition  changes,  such  as  phase  changes,  cause  variations  in 
the  amount  of  heat  needed  to  maintain  the  temperature  of  the  sample. 
These  variations  in  heat  energy  transfer,  or  heat  flux,  are  measured  in  a 
DSC  analysis  and  can  be  used  to  determine  material  composition  changes 
as  a  function  of  temperature.  This  test  was  also  performed  on  a  limestone 
sample  in  a  nitrogen  atmosphere  during  a  500°C  heat  treatment  with  a 
Netzch  Jupiter  Thermal  Analyzer.  Figure  3-7  shows  the  DSC  results. 


Figure  3-7.  DSC  results. 
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Negligible  changes  in  heat  flux  were  observed  that  indicated  no  phase 
changes  or  decomposition  at  temperatures  up  to  500C.  Significant  changes 
in  material  composition,  as  detected  by  DSC,  occur  at  temperatures  be¬ 
tween  700°C  and  iooo°C  and  can  cause  heat  flux  variations  of  approxi¬ 
mately  3  mW/mg.  These  results  along  with  those  from  DSC  and  TGA 
indicate  that  no  chemical  changes  are  happening  for  the  temperature 
range  of  interest.  Thus,  thermal  treating  is  a  viable  option  for  inducing  mi¬ 
crocracks  for  evaluation  with  mechanical  testing. 
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3.4  Quantifying  damage 

Several  methods  were  used  to  quantify  the  damage  from  the  heat  treat¬ 
ments  including  computed  tomography  (CT)  scanning,  scanning  electron 
microscopy  (SEM),  and  optical  microscopy  (OM).  Each  of  these  techniques 
will  be  discussed  in  this  section. 

3.4.1  Computed  tomography  scan 

CT  scanning  is  a  3-D  imaging  technique  that  uses  x-rays  to  produce  tomo¬ 
graphic  images  of  samples.  This  technique  was  performed  on  untreated, 
250°C,  and  450°C  heat  treated  samples  to  determine  if  it  can  be  used  to 
quantify  microstructural  changes  caused  by  heat  treating.  The  CT  scan  was 
performed  using  the  SkyScan  1173  high  energy  spiral  scan  micro-CT  with  a 
maximum  resolution  of  approximately  7  microns.  The  samples  used  in  this 
study  were  from  Frew  (2001)  and  were  cylinders  with  a  diameter  and 
length  of  12.7  mm  (0.5  in.)  These  smaller  samples  were  used  to  obtain  an 
image  resolution  of  approximately  8  microns;  larger  samples  could  only  be 
scanned  at  lower  resolutions. 

The  CT  scans  did  not  have  the  resolution  needed  to  capture  microcracking, 
so  the  images  were  used  to  determine  percent  void  increases  between  un¬ 
treated,  250°C,  and  450°C  heat  treated  samples.  Figure  3-8  shows  CT 
cross-section  images  for  each  sample  type.  For  the  percent  voids  calcula¬ 
tion,  the  void  area  was  determined  by  setting  the  grey  threshold  such  that 
the  undamaged  void  area  matched  the  porosity  of  the  virgin  material  (ap¬ 
proximately  14.0  percent).  The  same  threshold  was  used  for  the  250°C  and 
450°C  treated  samples  and  the  void  percentage  was  calculated.  The  CT  re¬ 
sults  show  the  250°C  and  450°C  treated  samples  had  void  area  increases 
to  16.8  percent  and  17.7  percent,  respectively. 


Figure  3-8.  CT  scan  cross-section  images  from  (Lto  R)  untreated,  250 °C,  and  450 °C 

heat  treated  samples. 
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The  increase  in  void  area  with  higher  temperature  heat  treatments  seen  in 
the  CT  scan  results  do  not  seem  reasonable  because  the  resolution  of  the 
CT  scans  is  larger  than  the  cracks,  which  should  be  the  major  contributor 
to  the  increase  in  void  area.  The  image  analysis  is  highly  sensitive  to  the 
chosen  threshold  value,  which  may  account  for  the  increases.  Until  CT 
scan  resolutions  of  1-2  microns  can  be  obtained  and  proper  threshold  val¬ 
ues  are  determined,  CT  scanning  is  not  recommended  for  use  when  deter¬ 
mining  crack  properties  and  void  areas  for  damage  quantification 
purposes. 

3.4.2  Scanning  Electron  Microscope 

SEM  was  used  in  the  initial  stages  of  the  microcrack  characterization  to 
provide  images  for  a  qualitative  analysis  of  the  changes  induced  by  the 
heat  treatments.  For  the  observations,  2-in.-diameter  samples  (same  as 
samples  used  in  sonic  velocity  tests)  were  heat  treated  to  500°C.  The  sur¬ 
faces  were  polished  with  diamond  paste  and  coated  with  gold-palladium. 
The  backscattered  SEM  images  were  taken  using  the  FEI  Nova  NanoSEM 
630.  The  SEM  images,  shown  in  Figure  3-9  and  Figure  3-10,  show  no  sig¬ 
nificantcracking  in  the  untreated  material.  Following  the  500°C  heat  treat¬ 
ment,  extensive  microcracking  along  grain  boundaries  and  inside  grains 
can  be  observed. 


Figure  3-9.  SEM  images  before  500  °C  heat  treatment. 
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Figure  3-10.  SEM  images  after  500  °C  heat  treatment. 


For  the  microstructural  characterization  in  this  work,  images  of  the  same 
locations  before  and  after  heat  treating  were  desired.  Though  SEM  images 
provide  sufficient  contrast  to  use  for  damage  quantification  purposes,  it 
was  unclear  how  the  gold-  palladium  coating  used  for  SEM  would  be  af¬ 
fected  by  the  heat  treatment  and  subsequent  characterization.  So,  optical 
microscopy  was  chosen  as  the  means  of  obtaining  images  for  microcrack 
quantification. 


3.4.3  Optical  microscope 

Optical  microscopy  was  performed  on  untreated  and  heat  treated  lime¬ 
stone  samples  to  provide  qualitative  and  quantitative  data  on  the  micro- 
structural  changes  as  a  function  of  the  maximum  temperatures  from  the 
heat  treatments.  OM  samples  were  created  by  coring  the  original  blocks 
with  a  conventional  rock  coring  barrel  and  were  cut  to  their  approximate 
final  length  with  a  slabbing  saw.  The  samples  were  then  polished  to  one 
micron  with  diamond  paste.  Samples  were  right  cylinders  with  a  length 
and  a  diameter  of  38.1  mm  (1.5  in.).  These  samples  were  slightly  smaller 
than  other  tests  samples  to  allow  them  to  fit  inside  the  polisher.  After  pol¬ 
ishing,  a  reference  mark  was  etched  in  each  sample  to  orient  the  micro¬ 
scope  in  the  same  locations  for  the  images  taken  before  and  after  heat 
treatments.  Only  surface  cracks  were  observed  because  crack  patterns 
formed  on  the  surface  show  only  minor  differences  from  interior  crack 
patterns  (Nolen-Hoeksema  and  Gordon  1987).  The  samples  were  exam¬ 
ined  using  a  Zeiss  Axiovert  200  optical  microscope  with  a  magnification 
range  of  50-2500  equipped  with  a  camera.  For  this  study,  images  were 
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taken  at  loox  magnification.  This  magnification  range  should  allow  cracks 
of  approximately  l  micron  in  width  to  be  observed. 

Figure  3-11  shows  the  limestone  before  and  after  250  °C  heat  treatment, 
and  Figure  3-12  shows  the  limestone  before  and  after  450  °C  heat  treat¬ 
ment.  The  figures  show  that  the  primary  effect  of  heat  treatment  in  this 
temperature  range  is  the  formation  of  new  microcracks.  The  samples 
heated  to  250  °C  introduced  cracks  into  the  sample  that  were  mostly  along 
grain  boundaries  with  few  intragranular  cracks.  The  450  °C  heat  treatment 
lengthened  some  existing  grain  boundary  cracks  and  introduced  many  in¬ 
tragranular  cracks  into  the  limestone  samples.  The  intragranular  cracks 
appear  to  develop  along  planes  of  crystallographic  weakness,  such  as 
cleavage  planes  and  pre-existing  flaws.  The  increasing  number  of  cracks 
with  increasing  temperature  results  in  the  decrease  in  the  measured  wave 
velocities. 

Figure  3-11.  Salem  limestone  before  (left)  and  after  (right)  250  °C  heat  treatment. 
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Quantification  of  the  induced  damage  with  OM  consisted  of  measuring  the 
microcrack  densities  for  the  250  °C  and  450  °C  heat  treatments.  The  quan¬ 
tification  was  performed  for  the  250  °C  and  450  °C  heat  treated  samples, 
because  these  are  the  heat  treatments  used  on  the  mechanical  test  samples 
discussed  in  the  next  chapter.  A  stereological  technique,  or  geometric 
probability  approach  (Fredrich  and  Wong  1986;  Wong  1985),  was  used  to 
quantify  the  microcrack  density.  The  approach  consists  of  counting  the 
number  of  crack  intersections  with  a  parallel  array  of  equally  spaced  lines. 
For  this  work,  an  array  of  10  parallel  lines  of  1  mm  length  spaced  0.1  mm 
apart  was  overlaid  on  the  OM  images.  Thus,  a  test  area  of  10  mm2  was  cov¬ 
ered  for  each  image.  A  total  of  128  images  were  analyzed  resulting  in  a  to¬ 
tal  test  area  of  1280  mm2.  The  OM  images  were  taken  in  two  locations  on 
the  sample,  the  center  and  midway  between  the  center  and  edge  of  the 
sample.  The  criterion  for  a  crack  was  that  it  has  to  be  visible  at  loox  mag¬ 
nification  and  that  it  has  an  aspect  ratio  greater  than  about  1:10  to  restrict 
the  inclusion  of  pores. 

The  stereological  data  for  Salem  limestone  heat  treated  to  250  °C  and 
450  °C  is  shown  in  Table  3-2,  where  the  stereological  parameter,  Pl,  is  the 
average  number  of  cracks  per  unit  length.  The  undamaged  material  had  a 
negligible  amount  of  cracks  (about  1  crack  per  10  mm2  test  array  area  that 
usually  did  not  intersect  a  test  array  line).  The  microcrack  density  in¬ 
creased  with  increasing  temperatures  with  the  450  °C  heat  treated  samples 
having  more  than  three  times  the  microcrack  density  than  the  250  °C  heat 
treated  samples. 
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Table  3-2.  Stereological  data  from  OM  images  for  heat  treated  Salem  limestone. 


Sample  Type 

Test  Area,  mm2 

Pl  Total, /mm 

Untreated 

640 

Negligible 

250°C  Heat  Treatment 

320 

0.94 

450°C  Heat  Treatment 

320 

3.14 

3.5  Summary 

This  chapter  described  the  limestone  material  used  in  this  work  and  the 
experimental  methods  used  to  induce  and  quantify  microstructural 
changes  in  the  material.  Thermal  treating  was  chosen  to  induce  mi¬ 
crocracks  in  the  material  rather  than  freezing/thawing  or  mechanical  test¬ 
ing  because  it  is  inexpensive  and  has  less  variables  to  consider.  Wave 
velocities  were  measured  as  a  non-destructive  means  of  evaluating  the  ef¬ 
fects  of  damage.  Wave  velocities  decreased  with  increasing  maximum  tem¬ 
peratures  used  in  the  heat  treatments  indicating  that  the  induced 
microstructural  changes  were  a  function  of  that  maximum  temperature. 
The  wave  velocities  also  showed  little  variability  for  a  particular  heat  treat¬ 
ment,  which  showed  that  the  heat  treatments  could  induce  consistent  lev¬ 
els  of  microcracking.  XRD,  DSC,  and  TGA  were  performed  to  provide  a 
range  of  temperatures  in  which  to  heat  the  samples  without  changing  the 
material  composition.  The  results  of  those  three  analyses  confirmed  that 
no  composition  changes  occurred  at  500  °C  or  below;  therefore,  room  tem¬ 
perature  to  500  °C  became  the  temperature  range  of  interest  for  the  mi¬ 
crocrack  characterization  and  mechanical  testing. 

CT  scanning,  SEM,  and  OM  were  used  to  observe  microstructural  changes 
caused  by  the  250  °C  and  450  °C  heat  treatments.  CT  scans  were  able  to 
provide  3D  images  of  the  untreated  and  treated  samples.  However,  the 
resolution  of  the  CT  scan  was  too  large,  -7  microns,  to  provide  microcrack 
information,  which  requires  a  resolution  of  around  1-2  microns.  An  SEM 
analysis  was  also  performed  to  give  an  initial  indication  of  microcracking 
caused  by  the  heat  treatments.  For  the  microcrack  characterization,  im¬ 
ages  in  the  same  locations  before  and  after  heat  treating  were  preferred. 
Because  SEM  requires  a  gold-palladium  coating  in  order  to  obtain  images, 
OM  analysis  was  chosen  for  the  microstructural  characterization  to  avoid 
any  effects  the  coating  may  have  on  the  sample  during  heating.  The  OM 
analysis  yielded  both  qualitative  and  quantitative  data  on  the  microstruc¬ 
tural  changes  associated  with  the  heat  treatments.  OM  images  showed  an 
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increase  in  grain  boundary  and  intragranular  cracks  with  increasing  maxi¬ 
mum  temperatures  in  the  heat  treatments.  A  stereological  technique  pro¬ 
vided  microcrack  densities  for  the  corresponding  heat  treatments.  The 
microcrack  density,  as  given  by  a  stereological  parameter,  increased  as  a 
function  of  the  maximum  temperature  in  the  heat  treatments. 
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4  Mechanical  Testing 

4.1  Overview 

Mechanical  tests  were  performed  to  characterize  the  mechanical  response 
of  the  intact  and  damaged  limestone  at  low  strain  rates  for  various  stress 
states  and  high  strain  rates  for  uniaxial  compression.  This  section  de¬ 
scribes  the  test  methods  and  results.  The  quasi-static  tests  using  servo-hy¬ 
draulic  loading  will  be  discussed  followed  by  the  dynamic  tests  using  the 
split  Hopkinson  pressure  bar  (SHPB).  In  this  chapter,  the  samples  heat 
treated  to  a  maximum  temperature  of  250°C  and  450°C  will  be  referred  to 
as  low  damage  and  high  damage  samples,  respectively.  Low  damage  and 
high  damage  are  also  considered  synonymous  with  low  and  high  levels  of 
microcracking. 

4.2  Quasi-static  testing 

Quasi-static  tests  were  performed  on  undamaged  and  damaged  samples  to 
determine  how  damage  affects  various  mechanical  properties.  A  total  of  45 
compression  tests  were  successfully  performed  including  9  unconfined 
compression  (UC)  tests,  24  triaxial  compression  (TXC)  tests,  3  hydrostatic 
compression  (HC)  tests,  and  9  uniaxial  strain/constant  volume  (UX/CV) 
strain  path  tests.  All  of  the  tests  were  performed  quasi-statically  with 
strain  rates  ranging  from  10  5  to  io_4  sA  At  these  rates,  times  to  reach  peak 
load  varied  from  5  to  30  minutes. 

Prior  to  performing  the  mechanical  tests,  the  height,  diameter,  and  weight 
of  each  sample  was  measured  and  used  to  compute  the  wet  density  of  ma¬ 
terial.  Following  the  test,  post-test  water  contents  were  measured  accord¬ 
ing  to  ASTM  D  2216  (ASTM  2005d).  Using  the  post-test  water  content, 
wet  density,  and  grain  density  (determined  to  be  2.75  Mg/m3),  porosity, 
degree  of  saturation,  and  air,  water,  and  solid  volumes  were  calculated. 
The  average  values  for  the  limestone  samples  for  wet  density  (which 
equals  the  dry  density  because  the  samples  were  dried  before  testing)  and 
porosity  were  2.34  Mg/m3  and  14.9  percent,  respectively.  Values  for  each 
sample  can  be  found  in  the  Appendix. 
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4.2.1  Sample  preparation 

Samples  for  the  quasi-static  mechanical  tests  were  cored  using  a  50.8  mm 
(2  in.)  diameter,  diamond-bit  core  barrel  following  the  procedures  pro¬ 
vided  in  ASTM  C  42  (ASTM  2005b).  The  samples  were  cut  to  the  correct 
length,  and  the  ends  were  ground  flat  and  parallel  to  each  other  (within 
±0.025  mm)  and  perpendicular  to  the  core  side  following  the  procedures 
provided  in  ASTM  D  4543  (ASTM  2005c).  The  finished  samples  had  a 
nominal  diameter  and  length  of  50.8  mm  (2  in.)  and  114.0  mm  (4.5  in.), 
respectively. 

Hardened  steel  caps  were  placed  on  the  ends  of  each  sample  before  test¬ 
ing. 


Except  for  the  unconfined  compression  tests,  two  0.6  mm  thick  latex 
membranes  were  placed  around  the  samples  (see  Figure  4-1),  and  the  exte¬ 
rior  of  the  outer  membrane  was  coated  with  liquid  synthetic  rubber  to  pre¬ 
vent  deterioration  from  the  confining  fluid,  a  mixture  of  kerosene  and 
hydraulic  oil.  After  the  membrane  and  coating  were  applied,  the  samples 
were  placed  on  the  testing  device  and  the  instrumentation  was  attached. 


Figure  4-1.  Typical  test  sample  setup. 
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4.2.2  Test  instrumentation 

A  vertical  deflection  measurement  system,  consisting  of  two  linear  variable 
differential  transformers  (LVDTs),  was  used  for  all  tests  to  measure  the  ax¬ 
ial  deformation  of  the  samples.  These  were  placed  between  the  top  and 
base  caps  of  the  sample.  For  confined  tests,  a  linear  potentiometer  was 
mounted  to  the  outside  of  the  pressure  vessel  to  measure  the  displacement 
of  the  piston  applying  the  axial  load.  This  provided  an  additional  measure¬ 
ment  of  the  axial  displacement  in  case  the  LVDTs  malfunctioned  or  ex¬ 
ceeded  their  calibrated  range. 

A  radial  deflection  measurement  system,  consisting  of  lateral  deformeters, 
was  used  to  measure  the  radial  deflection  in  the  samples.  On  the  sides  of 
each  sample,  two  small,  steel  footings  were  attached  (see  Figure  4-1).  Each 
footing  had  a  threaded  post  that  extended  through  the  latex  membrane  on 
which  steel  caps  were  screwed  to  seal  the  membranes  to  the  footing.  The 
lateral  deformeters  were  attached  to  these  steel  caps  with  set  screws  (see 
Figure  4-2). 


Figure  4-2.  Lateral  deformeters  mounted  on  test  specimen. 


Outputs  from  the  instrumentation  were  electronically  amplified  and  fil¬ 
tered,  and  the  signals  were  recorded  by  computer-controlled,  16-bit,  ana- 
log-to-digital  converters.  The  data  acquisition  system,  sampling  at  a  rate  of 
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0.2  to  l  Hz,  converted  the  voltage  to  engineering  units  (displacement,  load, 
and  pressure)  and  stored  the  data  for  further  processing. 

4.2.3  Stress  and  strain  definitions 

During  the  mechanical  tests,  the  axial  and  radial  deformations,  the  axial 
load,  and  the  confining  pressures  (except  for  UC)  were  recorded.  These 
measurements  and  the  pretest  height  and  diameter  of  the  samples  were 
used  to  calculate  the  true  stresses  and  engineering  strains.  Compressive 
stresses  and  strains  are  considered  positive  for  these  tests. 

Axial  strain,  ea,  was  calculated  by  dividing  the  axial  deformation,  Ah,  by 
the  original  height,  ho.  Radial  strain  was  calculated  by  dividing  the  radial 
deformation,  Ad,  by  the  original  diameter,  do.  Volumetric  strain,  ev,  was 
assumed  to  be  the  sum  of  the  axial  strain  and  two  times  the  radial  strain, 
as  given  in  Eq.  4.1. 


Sv  —  £a  +  2&r  (4.1) 

The  principal  stress  difference  (PSD),  q,  is  determined  by  dividing  the  ax¬ 
ial  load  by  the  cross-sectional  area  of  the  sample,  A,  which  is  equal  to  the 
original  cross-sectional  area,  A0,  multiplied  by  (1  -  £r)2  as  given  in  Eq.  4.2. 


C[  (<Ta  (Tr) 


Axial  Load 
A0(l-zr)2 


(4.2) 


Here,  oa  is  the  axial  stress  and  or  is  the  radial  stress.  For  uniaxial  compres¬ 
sion,  the  principal  stress  difference  is  simply  the  axial  stress.  The  mean 
normal  stress  (MNS),  p,  is  the  average  of  the  applied  principal  stresses, 
given  in  Eq.  4.3  for  cylindrical  samples. 


V  = 


2Cj- 

3 


(4.3) 


4.2.4  Unconfined  compression 

Nine  UC  tests  (3  undamaged,  3  low  damage,  and  3  high  damage)  were  per¬ 
formed  to  provide  shear  and  failure  data  on  the  undamaged  and  damaged 
limestone  samples.  For  the  unconfined  compression  tests,  the  axial  load 
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was  applied  using  a  3.3  MN  (750,000  lb-force)  loader.  The  load  applica¬ 
tion  was  manually  controlled.  Only  top  and  base  caps,  a  load  cell,  and  ver¬ 
tical  and  radial  deformeters  were  needed  for  these  tests. 


Figure  4-3  shows  the  typical  PSD-axial  strain  curves  for  the  unconfined 
compression  tests  performed  on  undamaged,  low  damage,  and  high  dam¬ 
age  samples  (see  appendix  for  all  mechanical  test  results).  Table  4-1  shows 
the  mean  peak  PSD’s  and  the  percentage  of  the  PSD’s  compared  to  the  un¬ 
damaged  materials.  At  the  beginning  of  the  unconfined  compression  tests 
on  the  treated  samples  (up  to  approximately  0.07  percent  axial  strain  for 
low  damage  and  0.23  percent  for  450  °C  high  damage  samples),  cracks 
were  closing  perpendicular  to  the  loading  axis  causing  an  apparent  ductil¬ 
ity.  That  apparent  ductility  correlated  to  an  increase  in  volume  caused  by 
the  induced  microcracks. 


Figure  4-3.  From  left  to  right,  principal  stress  difference  vs.  axial  strain  and  principal 
stress  difference  vs.  volumetric  strain  for  undamaged  and  damaged  samples  during 

uniaxial  compression. 
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Table  4-1.  Mean  peak  principal  stress  difference  during  unconfined  compression  for 

undamaged  and  damaged  samples. 


Sample  Type 

Mean  Peak  PSD,  MPa 

Percent  of  Undamaged  PSD 

Undamaged 

72.3 

100 

Low  Damage 

65.3 

90.3 

High  Damage 

45.7 

63.2 

Figure  4-3  also  shows  the  PSD-volumetric  strain  curves.  Initially  the  volu¬ 
metric  strains  were  positive,  indicating  compression.  The  mean  peak  volu¬ 
metric  strains  for  the  undamaged,  low  damage,  and  high  damage  samples 
were  0.08%,  0.12%,  and  0.27%,  respectively.  Once  the  peak  volumetric 
strains  were  reached,  negative  volumetric  strains  were  produced,  because 
the  samples  were  expanding  more  quickly  in  the  radial  direction  than  the 
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axial  direction.  This  radial  expansion  was  due  to  crack  coalescence  and  di¬ 
lation.  Because  the  heat  treated  samples  had  a  larger  crack  density  ini¬ 
tially,  the  cracks  coalesced  and  expanded  more  quickly,  which  led  to  the 
onset  of  dilatation  at  lower  stresses. 

4.2.5  Triaxial  compression 

Twenty-four  triaxial  compression  tests  were  performed  to  provide  addi¬ 
tional  shear  and  failure  data.  Six  tests  each  were  performed  at  10  MPa,  20 
MPa,  100  MPa,  and  400  MPa  confining  pressures.  For  each  set  of  six  tests 
at  each  confining  pressure,  two  were  performed  on  undamaged,  two  on 
low  damage,  and  two  on  high  damage  samples.  In  these  tests,  the  axial 
load  was  applied  using  an  8.9  MN  (2,000,000  lb-force)  loader,  and  the 
confining  pressures  were  applied  with  a  600  MPa  capacity  pressure  vessel 
(see  Figure  4-4.  The  application  of  load,  pressure,  and  axial  displacement 
were  controlled  by  a  servo-controlled  data  acquisition  system. 


Figure  4-4. 600-MPa  pressure  vessel  details. 


The  TXC  tests  were  performed  in  two  stages,  hydrostatic  loading  followed 
by  shear  loading.  In  the  initial  stage,  the  samples  are  subjected  to  a  hydro¬ 
static  loading  up  to  the  desired  pressure,  e.g.  10  MPa,  20  MPa,  100  MPa, 
and  400  MPa  for  this  work.  This  stage  of  loading  measures  the  pressure- 
volume  response  of  the  material,  which  can  be  used  to  calculate  the  mate¬ 
rial’s  bulk  modulus.  Following  the  hydrostatic  loading  stage,  shear  loading 
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is  applied  by  holding  the  confining  pressure  constant  while  increasing  the 
axial  load.  For  this  work,  the  peak  strength  of  the  material  is  defined  as  the 
maximum  PSD  that  a  sample  can  support  or  the  PSD  at  15  percent  axial 
strain  during  shear  (whichever  comes  first).  This  stage  of  loading  provides 
the  PSD-axial  strain  response,  which  can  be  used  to  calculate  the  elastic 
moduli. 


Figure  4-5  to  Figure  4-8  show  the  PSD-axial  strain  curves  for  the  TXC  tests 
at  the  various  confining  pressures.  Table  4.2  shows  the  mean  peak  PSD  for 
the  TXC  tests.  Note  that  the  axial  strains  are  given  minus  the  axial  strains 
induced  during  the  hydrostatic  loading.  These  results  show  that  for  the 
range  of  confining  pressures  tested  the  strengths  did  not  change  signifi¬ 
cantly.  This  shows  that  for  some  confining  pressure  less  than  10  MPa, 
there  exists  a  threshold  pressure  after  which  damage  does  not  lower  the 
strength  of  the  material. 


Figure  4-5.  (Lto  R)  Principal  stress  difference  vs.  axial  strain  and  principal  stress 
difference  vs.  volumetric  strain  for  undamaged  and  damaged  samples  during  triaxial 
compression  with  10  MPa  confining  pressure. 
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Figure  4-6.  (Lto  R)  Principal  stress  difference  vs.  axial  strain  and  principal  stress 
difference  vs.  volumetric  strain  for  undamaged  and  damaged  samples  during  triaxial 
compression  with  20  MPa  confining  pressure. 
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Figure  4-7.  (Lto  R)  Principal  stress  difference  vs.  axial  strain  and  principal  stress 
difference  vs.  volumetric  strain  for  undamaged  and  damaged  samples  during  triaxial 
compression  with  100  MPa  confining  pressure. 
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Figure  4-8.  (Lto  R)  Principal  stress  difference  vs.  axial  strain  and  principal  stress 
difference  vs.  volumetric  strain  for  undamaged  and  damaged  samples  during  triaxial 
compression  with  400  MPa  confining  pressure. 
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Figure  4-5  to  Figure  4-8  also  show  the  PSD-volumetric  strain  curves  for 
the  TXC  tests  at  the  various  confining  pressures.  Initially,  the  volumetric 
strains  were  compressive  and  remain  positive  until  dilation  occurs  in  the 
samples.  Table  4-2  shows  the  mean  peak  volumetric  strain  during  shear 
for  the  TXC  tests. 


Table  4-2.  Mean  peak  PSD  and  volumetric  strain  for  TXC  tests  on  undamaged  and 

damaged  samples. 


Sample  Type 

Mean  Peak  PSD,  MPa 

Mean  Peak  Volumetric  Strain,  Percent 

TXC  at  10  MPa 

Undamaged 

87.0 

0.122 

Low  Damage 

89.2 

0.128 

High  Damage 

85.5 

0.289 
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Sample  Type 

Mean  Peak  PSD,  MPa 

Mean  Peak  Volumetric  Strain,  Percent 

TXC  at  20  MPa 

Undamaged 

113.0 

0.128 

Low  Damage 

124.6 

0.162 

High  Damage 

106.8 

0.240 

TXC  at  100  MPa 

Undamaged 

240.8 

1.64 

Low  Damage 

234.3 

1.07 

High  Damage 

241.7 

1.43 

TXC  at  400  MPa 

Undamaged 

572.9 

3.32 

Low  Damage 

567.2 

2.63 

High  Damage 

570.3 

3.56 

4.2.6  Hydrostatic  compression 

Three  hydrostatic  compression  tests  were  performed  up  to  a  pressure  of 
400  MPa  to  provide  undrained  compressibility  data  for  the  limestone  sam¬ 
ples.  One  test  was  performed  each  on  the  undamaged,  low  damage,  and 
high  damage  samples.  The  test  setup  was  the  same  as  that  used  in  the  ini¬ 
tial  stage  of  the  TXC  tests. 

The  pressure-volume  response  of  the  undamaged  and  damaged  samples 
can  be  seen  in  Figure  4-9.  During  the  tests,  the  pressure  was  intentionally 
held  constant  for  a  brief  period  (approximately  4  min.)  before  unloading. 
During  this  time  the  strains  continued  to  increase,  which  indicated  that 
the  material  is  susceptible  to  creep  at  high  pressures.  For  each  material, 
yield  occurs  and  plastic  strains  begin  to  accumulate  at  a  pressure  of  ap¬ 
proximately  175  MPa.  The  high  damage  sample  initially  had  a  lower  bulk 
modulus  (-2000  MPa)  but  began  to  stiffen  with  increased  strain.  After 
that  initial  stiffening,  the  damaged  material  had  a  bulk  response  similar  to 
that  of  the  undamaged  material  (-14000  MPa). 
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Figure  4-9.  Pressure  volume  response  for  undamaged  and  damaged  samples  during 

hydrostatic  compression  to  400  MPa. 
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4.2.7  Uniaxial  strain/constant  volume  strain  path  test 

Nine  UX/CV  tests  were  performed  on  the  samples  to  obtain  failure  and 
one-  dimensional  compressibility  data.  The  axial  load  was  applied  using 
an  8.9  MN  (2,000,000  lb-force)  loader,  and  the  confining  pressures  were 
applied  with  a  600  MPa  capacity  pressure  vessel.  Three  tests  were  per¬ 
formed  (one  for  each  damage  level)  to  100  MPa  confining  pressure  during 
the  UX  portion  of  the  test,  three  to  150  MPa  confining  pressure,  and  three 
to  200  MPa  confining  pressure.  These  tests  were  performed  in  two  stages, 
uniaxial  strain  loading  and  loading  that  maintains  constant  sample  vol¬ 
ume.  The  UX  portion  of  the  test  was  conducted  by  applying  an  axial  load 
and  confining  pressure  simultaneously  so  that  the  diameter  of  the  sample 
remains  unchanged  as  the  axial  load  is  applied.  Following  the  UX  loading, 
a  constant  axial-to-radial  strain  ratio  (ARSR)  of  -2.0  was  applied.  The 
ARSR  of  -2.0  insures  a  constant  volumetric  strain  as  the  sample  is  loaded. 
UX/CV  tests  are  useful  for  confirming  results  obtained  from  the  TXC  be¬ 
cause  the  data  approximately  follows  the  recommended  failure  surface  of 
the  material,  as  seen  in  Figure  4-10. 
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Figure  4-10.  Comparison  of  UX/CV  and  TXC  tests  with  recommended  failure  surface. 
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Figure  4-11  shows  the  pressure-volume  response  and  PSD  vs.  MNS  for  the 
UX/CV  tests  to  150  MPa  confining  pressure.  The  pressure-volume  re¬ 
sponse  shows  that  high  damage  samples  accumulate  more  strain  before 
reaching  the  same  pressures  as  the  less  damaged  samples.  This  response  is 
due  to  crack  closure  with  increased  loading.  The  vertical  records  following 
the  UX  portion  indicate  the  CV  portion  of  the  test  where  volumetric  strain 
remains  constant.  The  PSD-MNS  response  follows  the  failure  surface  of 
the  material  and  is  slightly  lower  than  the  TXC  results.  The  PSD-MNS  re¬ 
sponse  also  shows  that  increasing  levels  of  damage  requires  higher  axial 
loading  to  increase  the  pressure  during  UX.  For  the  CV  portions  of  these 
tests,  there  is  not  a  significant  change  in  the  failure  surface  between  the 
undamaged  and  damaged  material. 


Figure  4-11.  (L  to  R)  Mean  Normal  Stress  vs.  volumetric  strain  and  principal  stress 
difference  vs.  mean  normal  stress  for  undamaged  and  damaged  samples  during 

UX/CV  tests. 
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4.3  Dynamic  testing 

Dynamic  testing  was  performed  using  a  modified  split  Hopkinson  pressure 
bar  (SHPB).  For  a  conventional  compression  SHPB  (see  Figure  4-12),  a 
sample  is  placed  between  two  bars,  an  incident  (input)  and  transmitted 
(output)  bar,  made  of  an  elastic  material,  e.g.,  a  high  strength  steel.  A  third 
bar,  the  striker  bar,  is  propelled  (typically  by  a  gas  gun)  into  the  end  of  the 
incident  bar  generating  a  longitudinal  compressive  incident  wave.  When 
the  wave  reaches  the  sample,  part  of  it  passes  through  the  sample  and  the 
remainder  is  reflected  back  into  the  incident  bar.  The  elastic  displace¬ 
ments  caused  by  the  compressive  waves,  measured  by  strain  gages,  are 
used  to  determine  the  stress-strain  response  of  the  sample.  Figure  4-13 
shows  the  compression  SHPB  at  the  Center  for  Advanced  Vehicular  Sys¬ 
tems  (CAYS)  at  Mississippi  State  University  (MSU)  used  for  this  work. 


Figure  4-12.  Conventional  SHPB  setup. 
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Figure  4-13.  Compression  SHPB  setup  at  MSU  CAVS. 
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For  brittle  materials  such  as  limestone,  a  nondispersive  ramp  pulse  is  re¬ 
quired  to  maintain  dynamic  stress  equilibrium  and  a  constant  strain  rate 
for  the  duration  of  the  SHPB  tests  needed  for  valid  results  (Frew  2001).  A 
modified  SHPB  technique  that  uses  a  thin  copper  pulse  shaper  is  used  in 
this  work  to  obtain  the  required  ramp  pulse  (see  Figure  4-14).  Figure  4-15 
shows  the  incident  wave  for  a  test  without  a  pulse  shaper  and  the  incident 
wave  obtained  in  this  work.  The  tests  with  the  pulse  shaper  provide  a 
nearly  linear  ramp  in  the  loading  up  to  100  ps  and  eliminate  the  high  fre¬ 
quency  oscillations  that  appear  in  a  conventional  SHPB  test. 


Figure  4-14.  Modified  SHPB  Setup. 
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Figure  4-15.  Incident  bar  stress  wave  without  (left)  and  with  (right)  a  pulse  shaper. 
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The  bars  used  in  this  work  were  made  of  high-strength  steel  and  had  diam¬ 
eters  of  12.8  mm.  The  striker,  incident,  and  transmitted  bars  had  lengths 
of  203.2  mm,  2438.4  mm,  and  1219.2  mm,  respectively.  The  strain  gage  on 
the  incident  bar  was  located  1219  mm  from  the  sample,  and  the  transmit¬ 
ted  bar  gage  was  381  mm  from  the  sample.  A  3.97  mm  diameter,  0.79  mm 
thick  annealed  copper,  C11000  (Lewis  1979  in  Frew  et  al.  2001),  disk  was 
placed  on  the  impact  surface  of  the  incident  bar  as  a  pulse  shaper.  The 
samples  used  in  the  work  were  prepared  by  Frew  (2001)  and  a  limited 
amount  was  available. 


Figure  4-16  shows  the  stresses  in  the  samples  at  the  incident  bar-sample 
interface  and  the  transmission  bar-sample  interface  for  the  undamaged, 
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low  damage,  and  high  damage  samples.  Because  the  incident  and  trans¬ 
mission  bar  interface  stresses  are  in  good  agreement,  the  samples  are  ap¬ 
proximately  in  a  state  of  dynamic  equilibrium.  If  the  interface  stresses 
were  not  in  agreement  (e.g.,  the  transmission  bar  stress  diverges  and  sig¬ 
nificantly  exceeds  the  incident  bar  stress),  that  would  indicate  that  the 
sample  was  not  in  stress  equilibrium  resulting  in  invalid  results.  For  the 
undamaged  sample,  the  strain  rate  is  relatively  constant  until  the  sample 
begins  to  fail  at  approximately  90  ps.  As  the  damage  increases,  a  greater 
variation  in  the  measured  strain  rate  can  be  seen.  There  is  also  a  larger 
variation  in  the  incident-  and  transmission-bar  interface  stresses  at  the  in¬ 
itial  portion  of  the  test  data.  This  most  likely  occurs  because  a  greater 
number  of  cracks  perpendicular  to  the  applied  stresses  are  being  closed  in 
the  samples  with  induced  damage.  The  results  from  all  of  the  SHPB  tests 
can  be  seen  in  the  appendix. 


Figure  4-16.  Interface  stresses  and  strain  rate  from  modified  SHPB  test  on 
undamaged  (top  left),  low  damage  (top  right),  and  high  damage  samples  (bottom 
left),  and  a  comparison  of  the  stress-strain  responses  for  each  damage  level  (bottom 

right). 
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The  undamaged  limestone  sample  had  a  failure  stress  of  about  80  MPa  at 
a  strain  rate  of  80  s_1,  which  is  lower  than  the  failure  stress  of  the  low  and 
high  damage  samples  of  about  90  MPa.  This  is  probably  due  to  strength 
variations  between  the  individual  samples.  Because  of  the  small  sample 
sizes  (0.5  inch  length  and  diameter),  density  changes  between  samples 
also  play  a  large  role.  This  strength  is  also  significantly  lower  than  the  dy¬ 
namic  strength  determined  by  Frew  (2001)  of  about  120  MPa  for  a  strain 
rate  between  100  and  120  s1.  This  could  be  explained  by  several  factors  in¬ 
cluding  differences  in  strain  rates  used  (80  s_1  vs.  100-120  s_1),  differences 
in  equipment  (Frew’s  tests  were  done  at  the  U.S.  Army  Engineering  Re¬ 
search  and  Development  Center),  or  simply  differences  between  samples 

The  damaged  and  the  undamaged  samples  had  failure  stresses  of  about  90 
MPa  for  strain  rates  between  100  and  130  s1.  Though  similar  incident 
stress  pulses  were  imparted  to  the  samples,  higher  strain  rates  were  ob¬ 
served  for  increasing  levels  of  damage.  This  is  most  likely  due  to  an  initial 
lower  stiffness  in  the  samples  (similar  to  that  seen  in  the  quasi-static  tests) 
during  crack  closure  offering  little  resistance  initially  that  allows  the  strain 
rate  to  increase  at  a  higher  rate.  This  same  behavior  could  also  be  respon¬ 
sible  for  undamaged  and  damage  samples  having  the  same  dynamic 
strengths. 

Though  the  overall  stiffness  of  the  damaged  samples  is  lower,  the  failure 
stresses  are  consistent  regardless  of  the  damage  levels,  as  seen  in  Figure 
4.16  (bottom  right).  More  dynamic  tests  are  needed  to  confirm  if  these 
trends  remain  for  a  larger  sample  size. 

4.4  Summary 

Quasi-static  and  dynamic  testing  was  performed  to  characterize  the  me¬ 
chanical  response  of  the  intact  and  damaged  limestone.  The  quasi-static 
tests  provided  results  for  tests  generating  several  stress  states  including 
uniaxial  compression,  triaxial  compression,  hydrostatic  compression,  and 
uniaxial  strain  /  constant  volume.  These  tests  provided  evidence  that  mi¬ 
crocracking  affected  limestone  by  lowering  peak  strengths  at  confining 
pressures  lower  than  10  MPa  and  lowering  the  stiffness  regardless  of  the 
stress  state  (though  greater  effects  can  be  seen  at  lower  pressures).  Mi¬ 
crocracking  did  not  affect  the  limestone’s  strength  at  pressures  greater 
than  10  MPa. 
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Modified  SHPB  tests  were  also  performed  to  determine  the  effects  of  mi¬ 
crocracking  on  the  dynamic  response  of  the  limestone.  These  tests  showed 
that  microcracks  lower  the  stiffness  of  the  material  even  at  higher  strain 
rates.  Microcracking  did  not  seem  to  have  an  effect  on  the  dynamic 
strength  of  the  limestone.  However,  a  limited  number  of  samples  were 
available  for  SHPB  tests,  and  more  dynamic  tests  are  needed  to  confirm  if 
these  trends  remain  for  a  larger  sample  size. 
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5  Material  Modeling 

5.1  Overview 

The  ability  to  predict  the  response  of  damaged  quasi-brittle  materials  has 
many  applications  in  civil  and  military  engineering  problems,  including 
penetration  and  blast  resistant  design  in  defense  structures  and  structural 
integrity  assessment  for  public  infrastructure.  In  addition  to  damage  ef¬ 
fects,  these  applications  commonly  involve  large  strains,  high  local  pres¬ 
sures,  and  high  strain  rates,  which  require  consideration  when  modeling 
these  types  of  events.  For  example,  Figure  5-1  shows  a  typical  projectile 
penetration  problem  where  pressures  around  the  tip  exceed  400  MPa  and 
large  areas  exceed  pressures  of  200  MPa. 


Figure  5-1.  Projectile  penetration  into  concrete  target. 


To  show  how  the  experimental  data  obtained  in  this  work  can  be  used  to 
improve  quasi-brittle  material  models,  the  Holmquist-Johnson-Cook 
(HJC)  model  (1993)  has  been  examined  and  modified  based  on  the  dam¬ 
aged  material  response  observed  in  the  previous  chapter.  The  HJC  model 
was  chosen  because  it  is  widely  used  for  blast  and  impact  problems  and  it 
is  a  relatively  simple  model  that  takes  into  account  the  most  important  is¬ 
sues  of  brittle  material  behavior,  such  as  pressure  dependency,  rate  de¬ 
pendency,  pressure-volume  changes,  and  damage.  Some  aspects  of  quasi- 
brittle  material  behavior  are  not  implemented  in  the  HJC  model  such  as 
stiffness  degradation  and  stress-state  dependency  and  have  been  ad¬ 
dressed  in  other  studies  (Riedel  et  al.  1999;  Gebbeken  and  Ruppert  2000; 
Polanco-Loria  et  al.  2008).  For  this  work,  modifications  were  made  to  ac¬ 
count  for  shear  modulus  degradation  and  failure  surface  changes,  which 
were  the  main  effects  of  microcracking  on  the  material  behavior  observed 
in  the  experimental  data. 
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5.2  Holmquist-Johnson-Cook  model 

The  HJC  model  is  an  elastic-viscoplastic  model  with  isotropic  damage  for 
concrete  subjected  to  large  strains  and  high  strain  rates  and  pressures.  In 
the  HJC  model,  the  deviatoric  response  of  the  material  is  determined  by 
the  constitutive  relation  given  in  Eqs.  5.1  and  5.2.  The  material  is  linear 
elastic  until  the  failure  surface,  or  shear  limit  surface,  is  reached.  It  in¬ 
cludes  a  scalar  damage  term  that  affects  the  material  by  reducing  the  cohe¬ 
sive  strength  value,  A,  in  the  failure  surface. 

a*  =  [A(l  -  D)  +  BP*n]  [1  +  C  Ins*]  for  a*  <  SMAX  (5.1) 

a*  =  SMAX  for  a*  >  SMAX  (5.2) 

In  Eq.  (1),  o*  and  P*  are  the  deviatoric  stress  and  pressure  normalized 
with  respect  to  the  compressive  strength  of  the  material,/^,  respectively, 
and  £*  is  the  strain  rate  normalized  to  a  reference  strain  rate  (usually  1  s_ 
J).  Using  normalized  pressure  and  tensile  hydrostatic  pressure  values  al¬ 
lows  for  the  same  values  to  be  used  for  different  unit  systems.  B,  N,  and  C 
are  material  constants,  where  B  and  N  affect  the  pressure  dependency, 
and  C  affects  the  strain  rate  dependency.  SMAX  is  the  maximum  strength 
and  acts  as  a  limiting  surface.  D  represents  damage. 

A  scalar  damage  formulation  is  used  in  the  HJC  model  and  is  defined  as 
the  accumulation  of  equivalent  plastic  strain  and  plastic  volumetric  strain, 
given  in  Eq.  5.3. 


D  = 


V,  (A£p+A/ip) 
1 


(5.3) 


In  Eq.  2,  Asp  and  Ajip  are  the  equivalent  plastic  strain  increment  and  equiv¬ 
alent  volumetric  strain  increment,  respectively,  and  sfp  and  jifp  those 
strains  at  failure.  The  plastic  strains  at  failure,  efp  and  //p,  are  expressed  as 

si  +  nfp  =  D1(P*  +  T*)°2  >  EFMIN  (5.4) 


where  Di  and  D2  are  constants  and  P*  and  T*  are  the  normalized  pressure 
and  normalized  tensile  hydrostatic  pressure,  respectively.  EFMIN  is  used 
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as  a  lower  limit  on  the  failure  strain  to  prevent  fracture  from  low  magni¬ 
tude  tensile  waves.  Figure  5-2  illustrates  the  failure  surface  used  in  the 
HJC  model.  The  effect  of  damage  in  this  formulation  is  to  translate  the 
failure  surface  down,  and  no  shape  change  is  possible. 

Figure  5-2.  The  HJC  failure  surface  showing  the  effects  of  damage. 


P* 


The  HJC  model  describes  the  pressure-volume  response  of  the  material 
with  three  distinct  regions.  In  the  first  region  (at  pressures  lower  than  the 
crush  pressure),  the  material  is  linear  elastic.  The  bulk  modulus  in  this  re¬ 
gion  is  given  by 


(5.5) 


where  Pcmsh  is  the  crush  pressure  and  pCmsh  is  the  volumetric  strain  at  the 
crush  pressure.  The  second  region  occurs  at  pressures  higher  than  the 
crush  pressure  and  lower  than  the  locking  pressure.  The  bulk  response  in 
this  region  is  determined  by  interpolating  the  responses  of  the  first  and 
third  region.  The  third  region  describes  the  relationship  for  fully  dense 
material,  where  all  of  the  air  voids  are  compressed  out  of  the  material.  In 
this  region,  pressure  is  a  non-linear  function  of  volumetric  strain,  given  in 
Eq.  5.6. 


P  =  K1/I+  K2/I2  +  K3/I3 


(5.6) 


In  Eq.  5.6,  the  modified  volumetric  strain,  p,  is  given  instead  of  the  volu¬ 
metric  strain  to  provide  a  relationship  to  the  volumetric  strain  at  locking, 
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(iiock  (the  beginning  of  the  third  region  in  the  pressure-volume  response). 
The  modified  volumetric  strain  is  given  in  Eq.  5.7. 


- _  f1~filock 


(5.7) 


Because  all  of  the  voids  are  gone,  large  pressure  increases  occur  for  very 
small  changes  in  volumetric  strain.  Figure  5-3  illustrates  the  pressure-vol¬ 
ume  response  used  in  the  HJC  model. 


Figure  5-3.  Pressure-Volume  Response  in  HJC. 


5.3  HJC  modifications 

The  experiments  performed  in  this  work  show  that  damage  affects  the 
strength  of  the  material  only  at  low  confining  pressures  (<  10  MPa).  As 
pressure  increases,  the  strength  of  the  damaged  material  increases  until  it 
is  equal  to  that  of  the  undamaged  material  (see  Table  4.2).  It  also  shows 
the  elastic  moduli  are  affected  by  damage  at  confining  pressures  less  than 
100  MPa.  Modifications  were  made  to  the  HJC  model  to  incorporate  these 
damage  effects  on  the  failure  surface  and  shear  modulus,  which  provides  a 
more  realistic  material  response  for  various  levels  of  damage.  The  modifi¬ 
cations  to  the  HJC  failure  surface  will  be  discussed  next  followed  by  the 
changes  to  the  shear  modulus  formulation. 
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5.3.1  Failure  surface 

In  the  first  term  of  the  HJC  failure  surface  (Eq.  5.1),  damage  reduces  the 
cohesive  strength,  A,  which  shifts  the  failure  surface  of  the  material  to  a 
damaged  state  as  shown  in  Figure  5.2.  Once  the  material  has  accumulated 
any  amount  of  damage,  it  cannot  return  to  the  undamaged  failure  surface. 
The  experiments  performed  in  this  work  indicate  that  damage  affects  the 
shear  strength  of  the  material  in  a  pressure  dependent  manner.  At  low 
confining  pressures  (less  than  10  MPa),  the  effects  are  significant;  as  pres¬ 
sure  increases,  damage  effects  on  the  shear  strength  diminish.  To  incorpo¬ 
rate  this  pressure  dependency,  a  modification  was  made  to  the  yield 
surface  formulation  that  allows  damage  to  reduce  the  strength  of  the  mate¬ 
rial  only  at  low  pressures.  Eqs.  5.8  and  5.9  show  the  new  failure  surface 
formulation,  which  resembles  the  original  formulation  with  an  additional 
multiplier  on  the  damage  constant. 


a*  =  [i4(l  —  D*)  +  BP*W][1  +  C  In e*]  for  a*  <  SMAX 


(5.8) 


(5.9) 


D*  is  the  new  pressure  dependent  damage  term,  P*rhreshoid  is  the  pressure 
after  which  damage  does  not  affect  the  shear  strength,  and  DC  is  a  damage 
constant  that  determines  how  quickly  the  effect  of  damage  decreases  with 
pressure.  Figure  5-4  illustrates  the  failure  surface  of  the  original  HJC  and 
the  failure  surface  proposed  in  this  work. 


Figure  5-4.  The  original  HJC  failure  surface  (left) 
and  proposed  HJC  failure  surface  (right). 
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5.3.2  Shear  modulus 

The  experimental  results  indicate  that  damage  affects  the  elastic  moduli  of 
the  material,  and  that  behavior  is  not  captured  by  the  original  HJC  formu¬ 
lation.  As  shown  in  the  hydrostatic  compression  tests  (discussed  in  Chap¬ 
ter  4),  the  initial  bulk  modulus  of  the  damaged  material  was  lower  than  the 
undamaged  material,  but  the  bulk  modulus  increased  with  increasing 
pressure  to  match  the  bulk  response  of  the  undamaged  material. 

Therefore,  no  damage  effects  on  the  bulk  moduli  are  considered  here. 
However,  damage  effects  on  the  shear  modulus  were  persistent  to  high 
pressures  so  those  affects  were  included  in  the  modifications  to  the  HJC 
model  in  this  work. 

To  incorporate  the  effects  of  damage  in  the  shear  modulus  formulation,  a 
linear  interpolation  scheme  (Gebbeken  and  Ruppert  2000;  Polanco-Loria 
et  al.  2008)  shown  in  Eq.  5.10  was  used. 

G  =  G0(1-D*)  +  GdD*  (5.10) 

G,  G0,  and  GD  represent  the  current,  initial,  and  completely  damaged  shear 
modulus,  respectively.  D*  was  used  instead  of  D  because  the  experimental 
results  also  showed  a  pressure  dependence  on  the  shear  modulus. 

5.4  Determination  of  model  parameters 

The  modified  HJC  model  was  implemented  in  the  solid  mechanics  code  Si¬ 
erra  (SAND2011-7597).  The  Sierra  implementation  of  the  original  HJC 
model  requires  the  identification  of  20  parameters:  the  initial  density;  two 
elastic  constants  (bulk  and  shear  modulus);  six  strength  constants  includ¬ 
ing  the  unconfined  compressive  and  tensile  strength,  cohesive  strength  A, 
pressure  hardening  coefficient  B,  pressure  hardening  exponent  N,  and 
maximum  shear  strength  SMAX;  one  rate  sensitivity  parameter  C;  three 
damage  constants  (Di,  D2,  and  EFMIN);  and  seven  constants  to  describe 
the  pressure-  volume  response  including  those  defining  the  elastic  region 
(Pcrush  and  Ucrush),  the  crushing  region  (Plock  and  Ulock),  and  the  lock¬ 
ing  region  (Ki,  K2,  and  K3).  The  modifications  to  the  HJC  model  require 
three  additional  parameters:  the  normalized  threshold  pressure  P* Threshold, 
another  damage  constant  DC,  and  a  fully  damaged  shear  modulus  Gd. 
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Simulations  were  performed  to  compare  with  the  undamaged  and  dam¬ 
aged  material  data.  For  the  undamaged  material,  the  initial  damage  was 
set  to  zero.  For  the  highly  damaged  material  (heated  to  450°C),  the  initial 
damage  was  set  to  0.65.  This  value  was  determined  using  the  damage  pa¬ 
rameter  proposed  by  MacKenzie  [25]  because  it  relates  damage  to  the 
shear  and  bulk  modulus  of  the  material.  The  MacKenzie  damage  parame¬ 
ter,  xpD,  is  shown  in  Eq.  (5.11). 


V D  = 


4  G(K-K) 
K(4G+3K) 


(5.11) 


The  parameters  K  and  G  in  Eq.  5.11  are  the  elastic  bulk  modulus  and  shear 
modulus  of  the  undamaged  material,  respectively.  The  parameter  K  is  the 
bulk  modulus  of  the  damaged  material.  For  all  problems,  the  damage  pa¬ 
rameter  is  between  o  and  1,  or  o  <.  ipD  <  1.  When  the  material  is  undam¬ 
aged,  K  =  K,  ijjd  =  o,  and  for  completely  damage  materials,  K  =  o,  ipo  =  1. 
For  this  work,  the  bulk  modulus  and  shear  modulus  were  determined  from 
wave  velocity  measurements  on  undamaged  and  damaged  samples.  Figure 
5-5  shows  the  MacKenzie  parameter  for  the  various  heat  treatments. 


Figure  5-5.  MacKenzie  damage  parameter  for  various  heat  treatments. 
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Fpr  the  strength  parameters,  data  from  the  unconfined  and  confined  com¬ 
pression  tests  were  used  to  find  the  find  the  values  of  A,  B,  and  N.  Because 
no  tension  tests  were  performed,  the  tensile  strength  was  assumed  to  be 
related  to  the  compressive  strength,  i.e.,/Y  =  6.7 /fc  (ACI  318,  2011). 
With  no  initial  damage,  the  original  and  modified  HJC  results  are  identical 
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(see  Figure  5-6).  With  high  initial  damage  (Do=o.65),  the  original  HJC 
failure  surface  is  shifted  down  because  of  the  reduction  in  the  cohesive 
strength.  The  modified  HJC  matches  the  failure  surface  at  both  lower  pres¬ 
sure,  where  the  strength  is  reduced,  and  higher  pressures.  SMAX  was 
simply  set  to  a  high  number  because  the  experiments  indicate  that  the  ma¬ 
terial  will  continue  to  pressure  harden  for  the  range  of  pressures  tested. 

Figure  5-6.  (left)  Model  fit  using  HJC  and  (right)  modified  HJC  for  undamaged  material 
(D0=0)  and  (b)  damaged  material  (D0=0.65). 
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The  three  damage  constants  (Di,  D2,  and  EFMIN  in  Eq.  5.4)  were  chosen 
in  a  manner  similar  to  that  of  Holmquist  et  al.  (1993),  which  used  cyclic 
uniaxial  compression  tests  to  determine  EFMIN.  Instead  of  cyclic  tests,  the 
unconfined  compression  tests  on  undamaged  and  damaged  samples  were 
used  to  define  an  assumed  failure  surface  seen  in  Figure  5-7.  This  method 
indicates  a  total  loss  of  strength  at  ePf=o.on,  so  EFMIN  was  set  equal  to 
0.011.  D2  was  chosen  to  equal  1.0,  which  assumes  the  plastic  fracture 
strain  increases  linearly  with  pressure.  Di  was  then  calculated  using  Eq. 

5.4  and  the  uniaxial  compression  data. 

Hydrostatic  compression  data  was  used  to  determine  the  constants  (Eqs. 

5.5  and  5.6)  that  define  the  pressure-volume  response  of  the  material.  No 
damage  effects  on  the  bulk  modulus  were  incorporated  into  the  modified 
version  of  the  HJC  so  the  pressure-  volume  response  is  the  same  for  the 
original  and  modified  HJC.  The  locking  region  in  the  pressure-volume  re¬ 
sponse  was  assumed  linear  for  the  purposed  of  this  work. 


ERDC/GSL  TR-17-15 


57 


Figure  5-7.  Uniaxial  compression  data  used  to  determine  HJC  damage  constants. 


After  fitting  all  of  the  constants  for  the  original  HJC  model,  the  three  addi¬ 
tional  constants  for  the  HJC  modifications  were  determined  including  the 
normalized  threshold  pressure  P* Threshold,  the  damage  constant  DC,  and  a 
fully  damaged  shear  modulus  Gd.  P* Threshold  was  set  to  the  brittle-to-ductile 
transition  of  the  material  to  allow  softening  of  the  material  to  occur  up  to 
the  brittle-to-ductile  transition.  The  damage  constant  DC  determines  how 
quickly  the  effects  of  damage  disappear  with  pressure  and  was  fit  to  give 
the  desired  results.  Lastly,  the  fully  damaged  shear  modulus,  Gd,  was  de¬ 
termined  from  the  sonic  velocities  of  the  material  heat  treated  to  the  high¬ 
est  temperature  (8oo°C). 

To  determine  the  strain  rate  parameter  C  (Eq.  5.1),  a  numerical  analysis  of 
the  SHPB  experiments  discussed  in  the  previous  chapter  was  performed. 
Initial  simulations  were  performed  with  no  strain  rate  dependence,  C=o. 
Then,  the  strain  rate  dependence  was  gradually  increased  until  the  simu¬ 
lated  strain-time  response  matched  that  of  the  experiment  (see  Figure 
5-8).  These  simulations  will  be  discussed  in  more  detail  later. 

Table  5-1  shows  the  rate-independent  and  dependent  constants  for  the 
modified  HJC  model  for  Salem  limestone. 
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Figure  5-8.  Strain-time  response  for  SHPB  simulations  with  and  without  rate 

dependence. 
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Table  5-1.  Constants  for  the  modified  HJC  model. 


Original  HJC  Constants 

No  Strain  Rate  Dep. 

Strain  Rate  Dep. 

Initial  Density,  g/cc 

2.30 

2.30 

Initial  Shear  Modulus,  MPa 

10000.00 

10000.00 

Cohesive  Strength  A 

0.50 

0.50 

Pressure  Hardening  Coefficient  B 

1.286 

1.373 

Pressure  Hardening  Exponent  N 

0.8362 

0.8273 

Strain  Rate  Coefficient  C 

0.00 

0.0038 

Compressive  Strength  F'c,  MPa 

72.00 

72.00 

Max  Stress 

13.00 

13.00 

Pcrush,  MPa 

195.00 

195.00 

Ucrush 

0.013 

0.013 

Pressure  Coefficient  Kl,  MPa 

44300.00 

44300.00 

Pressure  Coefficient  K2,  MPa 

0.00 

0.00 

Pressure  Coefficient  K3,  MPa 

0.00 

0.00 

Plocki,  MPa 

665.00 

665.00 

Ulock 

0.13 

0.13 

Max  Tensile  Pressure  T,  MPa 

6.00 

6.00 

Damage  Coefficient  D1 

0.05 

0.05 

Damage  Exponent  D2 

1.00 

1.00 

Minimum  Failure  Strain 

0.011 

0.011 

New  Constants 

P*threshold 

3.50 

3.50 

DC 

8.00 

8.00 

Gd,  MPa 

2000.00 

2000.00 
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5.5  Results 

Numerical  simulations  using  the  material  parameters  in  Table  5.1  were 
performed  with  the  original  HJC  and  modified  HJC  models  to  compare 
with  the  experimental  stress-  strain  data.  The  solid  mechanics  code  Sierra 
(SAND2011-7597)  was  used  for  both  the  quasi-static  and  dynamic  simula¬ 
tions. 

5.5.1  Quasi-static  test  analysis 

The  quasi-static  simulations  consisted  of  unconfined  compression  (UC) 
and  triaxial  compression  (TXC)  tests  and  were  analyzed  using  a  single 
brick  element.  The  results  labeled  undamaged  indicate  no  induced  damage 
in  the  experimental  samples,  and  Do=o  in  the  material  model  setup.  The 
results  labeled  high  damage  indicate  experimental  samples  with  damage 
induced  from  the  450°C  heat  treatment,  and  Do=o.6s  per  the  MacKenzie 
damage  parameter  mentioned  earlier.  Post-peak  data  is  given  for  compari¬ 
son  purposes  only;  the  experimental  and  numerical  implications  associ¬ 
ated  with  post-peak  response  are  not  discussed  because  it  is  beyond  the 
scope  of  this  analysis. 

Figure  5-9  compares  the  experimental  results  with  the  original  and  modi¬ 
fied  HJC  models  for  the  uniaxial  compression  tests.  The  principal  stress 
difference  (PSD)-axial  strain  response  is  given  for  the  comparisons.  The 
modified  HJC  provides  better  results  for  the  strength  and  stiffness  under 
uniaxial  compression.  The  lower  initial  stiffness  (due  to  crack  closure)  ob¬ 
served  in  the  experimental  results  (circled  in  Figure  5-9)  is  not  captured  by 
the  modified  HJC,  which  lead  to  differences  in  the  amount  of  axial  strain 
induced  by  the  loading.  Beyond  the  lower  initial  stiffness  in  the  experi¬ 
mental  results,  the  simulation  and  experimental  responses  are  approxi¬ 
mately  parallel. 
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Figure  5-9.  Unconfined  compression  data  compared  with  original  HJC  (left)  and 

modified  HJC  model  (right). 


0.0015 


0.003  0.0045  0.006 

Axial  Strain 


0.0075  0.009 


Figure  5-10  to  Figure  5-12  compare  the  experimental  results  with  the  origi¬ 
nal  and  modified  HJC  models  for  the  triaxial  compression  tests.  The  HJC 
results  show  that  a  reduction  of  the  failure  surface  with  increasing  damage 
does  not  represent  the  actual  effects  of  damage  especially  at  higher  pres¬ 
sures.  The  modifications  to  the  HJC  implemented  in  this  work  better  cap¬ 
ture  the  effects  of  damage  on  the  yield  surface,  which  provides  a  better 
prediction  of  the  peak  stresses  obtained  during  triaxial  loading.  The  modi¬ 
fied  HJC  also  captures  some  of  the  damage  effects  on  the  elastic  moduli. 


Figure  5-10.  Triaxial  compression  data  (10  MPa  confining  pressure)  compared  with 
original  HJC  (left)  and  modified  HJC  model  (right). 


Axial  Strain 


ERDC/GSL  TR-17-15 


61 


Figure  5-11  Triaxial  compression  data  (20  MPa  confining  pressure)  compared  with 
original  HJC  (left)  and  modified  HJC  model  (right). 
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Figure  5-12.  Triaxial  compression  data  (100  MPa  confining  pressure)  compared  with 
original  HJC  (left)  and  modified  HJC  model  (right). 
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5.5.2  Dynamic  test  analysis 

The  numerical  analysis  of  the  dynamic  tests  consisted  of  simulating  the 
split  Hopkinson  pressure  bar  (SHPB)  tests.  The  SHPB  geometry,  given  in 
Chapter  4,  was  modeled  in  quarter  symmetry  with  168,960  brick  elements, 
as  shown  in  Figure  5-13.  Instead  of  explicitly  modeling  the  striker  bar  and 
pulse  shaper,  the  stress  wave  created  by  the  striker  bar  in  the  experiments 
was  input  into  the  incident  bar.  The  incident  strain  measurements  from 
the  analysis  match  the  incident  bar  gage  measurements  from  the  experi¬ 
ments  indicating  that  explicitly  modeling  the  striker  bar  and  pulse  shaper 
is  not  necessary. 
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Figure  5-13.  SHPB  geometry  used  in  numerical  analysis. 


Figure  5-14.  SHPB  simulation  of  undamaged  and  damaged  samples  with  (left) 
original  HJC  and  (right)  modified  HJC  model. 
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Figure  5-14  compares  the  experimental  results  of  the  SHPB  tests  with  the 
original  and  modified  HJC  models.  For  the  initially  undamaged  samples, 
the  two  models  give  similar  responses.  For  the  damaged  samples,  the  orig¬ 
inal  HJC  (Figure  5-14  left)  over  predicts  the  elastic  stiffness  and  shows 
much  lower  strength,  approximately  57  percent  lower,  than  the  experi¬ 
mental  results.  The  modified  HJC  (Figure  5-14  right)  matches  the  yield 
strength  and  stiffness  of  the  experiments  better  than  the  original  HJC 
though  the  yield  strength  is  still  lower  than  the  experimental  results.  The 
modified  HJC  also  captures  the  higher  initial  stiffness  observed  in  the 
damaged  experimental  results  (seen  at  stresses  lower  than  10  MPa).  This 
may  indicate  that  the  higher  initial  stiffness  is  a  structural  response  in¬ 
stead  of  a  material  response. 

Furthermore,  Figure  5-14  illustrates  the  magnitude  of  the  under  prediction 
of  strain  energy  (area  under  the  stress-strain  curve)  by  the  HJC  model  for 
high  strain  rate  events  on  initially  damaged  materials.  Table  5-2  compares 
the  energies  from  the  experiments  and  two  models.  Though  the  models 
give  similar  results  for  the  undamaged  material,  the  HJC  model  yields  a 
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94  percent  decrease  from  the  strain  energy  in  the  experiment  on  the  ini¬ 
tially  damaged  material.  The  modified  HJC  performs  significantly  better 
with  a  53  percent  decrease. 


Table  5-2.  Energy  comparison  for  SHPB  tests  with  HJC  and  modified  HJC  models. 


Strain  Energy  (MPa) 

Damage  Level 

Experiment 

HJC 

Mod  HJC 

No  Damage 

0.427 

0.205 

0.213 

High  Damage 

0.644 

0.040 

0.301 

The  numerical  analysis  of  the  SHPB  tests  was  also  used  to  examine 
whether  the  strength  increase  observed  in  the  tests  was  due  solely  to  strain 
rate  or  if  radial  confinement,  caused  by  inertia  or  friction,  caused  the 
strength  to  increase.  During  a  SHPB  test,  the  rapid  axial  strain  accelera¬ 
tion  causes  radial  inertial  effects,  which  can  lead  to  radial  confining  pres¬ 
sures.  Friction  between  the  bar-sample  interfaces  also  causes  radial 
confinement  near  the  ends  of  the  sample.  In  pressure  dependent  materi¬ 
als,  the  confining  pressures  due  to  inertia  and  friction  could  cause  strength 
increases  that  may  be  misinterpreted  as  a  strain  rate  effect. 

To  investigate  the  inertia  and  friction  effects  during  the  SHPB  tests,  the 
strain  rate  independent  constants  (Table  5.1)  were  used  to  model  the 
SHPB  geometry  used  in  the  experiments.  Strain  rate  independent  models 
have  been  used  in  other  studies  to  investigate  the  effects  of  non-strain-rate 
effects  (Bertholf  and  Karnes  1975;  Meng  and  Li  2003;  Li  and  Meng  2003; 
Li  et  al.  2009)  because  any  differences  between  the  quasi-static  uncon¬ 
fined  behavior  and  the  SHPB  behavior  using  strain  rate  independent  pa¬ 
rameters  can  be  attributed  to  effects  unrelated  to  strain  rate. 

The  effect  of  inertia  on  the  sample  was  determined  by  simulating  the 
SHPB  test  with  no  friction  between  the  bar-sample  interfaces  and  using 
the  strain  rate  independent  constants.  Thus,  any  differences  between  the 
unconfined  compressive  strength  and  strengths  determined  from  these 
simulations  would  be  attributed  to  inertia  causing  radial  confinement.  Ta¬ 
ble  5-3  shows  the  results  of  the  simulations  on  the  initially  undamaged 
and  damaged  material.  The  insignificant  changes  between  the  UC  and 
SHPB  test  results  provide  evidence  that  radial  confinement  due  to  inertia 
did  not  increase  the  strength  of  the  sample  during  the  SHPB  test. 


ERDC/GSL  TR-17-15 


64 


Table  5-3.  Comparison  of  strengths  from  UC  simulations  and  SHPB  simulations  with 

no  friction  and  no  rate  dependency. 


Test  Type 

Peak  Sample 
Strength,  MPa 

Do=0 

Do=0.65 

UC 

74.28 

48.82 

SHPB 

74.55 

48.34 

To  illustrate  the  effects  of  inertia,  pressure  was  plotted  at  three  different 
times  during  the  SHPB  test.  Figure  5-15  shows  the  stress  levels  at  which 
these  plot  were  made.  Figure  5.16  (top)  shows  Point  A  in  Figure  5-15,  Fig¬ 
ure  5.16  (middle)  shows  Point  B,  and  Figure  5.16  (bottom)  shows  Point  C. 


Figure  5-15.  Stress-strain  response  from  SHPB  test  showing  stress  levels  for  images 

in  Figures  5.16  and  5.17. 
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Figure  5-16.  Pressure  in  sample  during  SHPB  test. 
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Because  inertia  did  not  have  an  effect  on  these  tests,  the  pressures  in  Fig¬ 
ure  5.16  are  simply  one  third  of  the  axial  stress  in  the  sample.  To  view  the 
effects  of  inertia  more  clearly,  one  third  of  the  axial  stress  was  subtracted 
from  the  pressure  (shown  in  Figure  5-17),  which  is  the  effective  confining 
pressure  throughout  the  SHPB  test.  The  effective  confining  pressure  re¬ 
mains  zero  throughout  the  entire  test  showing  no  additional  confining 
pressure  due  to  inertia.  Because  the  strain  rate  found  in  these  simulations 
is  approximately  80  s-1,  these  results  are  supported  by  the  findings  of  other 
studies  (Malvern  et  al.  1985;  Tang  et  al.  1992)  where  inertia  did  not  play  a 
significant  role  at  strain  rates  below  100  s1. 
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Figure  5-17.  Effective  confining  pressure  during  SHPB  test  showing 
inertial  effects  on  sample. 


To  determine  the  effects  of  friction  on  the  dynamic  strength,  a  parametric 
analysis  using  the  SHPB  simulations  was  performed  by  varying  the  con¬ 
stant  friction  coefficients  for  the  bar-sample  interfaces  from  o.o  to  0.5. 
This  range  was  chosen  because  it  is  the  range  of  friction  values  most  likely 
found  in  a  SHPB  test  (Bertholf  and  Karnes  1975).  In  this  analysis,  the 
strain  rate  independent  constants  were  also  used.  Table  5-4  shows  the  re¬ 
sults  of  the  analysis.  An  increase  of  two  percent  in  peak  sample  strength  is 
observed  from  no  friction  to  a  friction  coefficient  of  0.1  for  the  initially  un¬ 
damaged  material.  Peak  sample  strengths  increase  a  negligible  amount 
thereafter.  For  the  initially  damaged  sample,  the  strength  increases 
slightly  and  then  decreases  as  the  friction  coefficient  increases.  This  was 
due  to  localized  element  failure  that  happened  early  in  the  test  and  gov¬ 
erned  the  solution  despite  mesh  refinement.  Regardless,  negligible 
changes  in  the  strength  were  seen  with  varying  levels  of  friction;  however, 
further  discussion  will  be  limited  to  the  initially  undamaged  material  be¬ 
cause  of  this  numerical  issue. 
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Table  5-4.  Variation  in  peak  strength  with  various  constant  friction  coefficients 
between  bar-sample  interfaces  in  SHPB  simulation. 


Constant  Friction 
Coefficient 

Peak  Sample 
Strength,  MPa 

Do=0 

Do=0.65 

0.0 

74.55 

48.34 

0.1 

76.21 

49.55 

0.2 

76.37 

48.68 

0.3 

76.37 

47.47 

0.4 

76.37 

46.59 

0.5 

76.37 

46.96 

Similar  to  the  plots  made  for  the  inertia  analysis,  Figure  5.18  shows  pres¬ 
sure  at  different  times  during  the  SHPB  test  for  a  constant  pressure  coeffi¬ 
cient  of  0.1.  Figure  5-18  (top)  shows  Point  A  in  Figure  5-15,  Figure  5-18 
(middle)  shows  Point  B,  and  Figure  5-18  (bottom)  shows  Point  C.  In  con¬ 
trast  to  Figure  5.16  (where  friction  is  zero),  Figure  5-18  shows  a  pressure 
gradient  that  increases  as  the  SHPB  test  progresses. 


Figure  5-18.  Pressure  in  the  sample  during  the  SHPB  with  a  constant  friction 

coefficient  of  0.1. 


ERDC/GSL  TR-17-15 


68 


Figure  5-19  shows  corresponding  plots  of  pressure  minus  the  pressure  due 
to  axial  load  to  illustrate  the  effective  confining  pressures  during  the  test. 
Figure  5-19  (bottom),  taken  at  peak  axial  load,  indicates  that  these  confin¬ 
ing  pressures  can  reach  approximately  10  MPa  towards  the  center  of  the 
sample  and  5  MPa  along  the  outer  radius.  From  the  quasi-static  TXC  test 
with  10  MPa  confining  pressure,  an  increase  in  strength  of  approximately 
40  MPa  can  be  expected  in  the  areas  with  10  MPa  confinement.  These  con¬ 
fining  pressures  between  5  and  10  MPa  on  the  end  of  the  samples  resulted 
in  the  two  percent  increase  in  the  sample  strength. 

Figure  5-19.  Effective  confining  pressure  during  SHPB  test  with  a  constant  friction 
coefficient  of  0.1  showing  the  effects  of  friction. 


Though  the  confining  pressures  due  to  friction  in  these  tests  did  not  result 
in  a  significant  increase  in  strength,  it  is  an  indication  that  care  should  be 
taken  when  using  shorter  specimens,  which  would  place  a  larger  percent¬ 
age  of  the  sample  volume  under  the  increased  confining  pressures.  The 
presence  of  these  pressures  should  also  be  taken  intoconsideration  when 
interpreting  the  results  of  the  SHPB  tests  for  samples  with  microcracks. 
The  quasi-static  mechanical  tests  indicated  that  confining  pressures  as  low 
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as  10  MPa  were  enough  to  eliminate  the  effects  of  microcracking  on  the 
sample  strength.  Figure  5-19  shows  that  a  friction  coefficient  as  low  as  0.1 
is  enough  to  induce  pressures  that  would  alter  the  effects  of  microcracks 
on  a  portion  of  the  sample  during  SHPB  tests.  Thus,  this  provides  further 
justification  that  friction  should  be  eliminated  during  high  rate  testing. 
Friction  and  inertial  confinement  should  be  investigated  further  to  deter¬ 
mine  its  effects  on  microstructural  mechanisms  during  high  rate  loading. 

5.6  Summary 

The  HJC  model  has  been  examined  and  modified  based  on  the  damaged 
material  response  observed  in  the  previous  chapter.  Modifications  were 
made  to  account  for  shear  modulus  degradation  and  failure  surface 
changes,  which  were  the  main  effects  of  damage  on  the  material  behavior 
observed  in  the  experimental  data.  The  original  and  modified  HJC  models 
were  used  in  a  numerical  analysis  of  the  quasi-static  and  dynamic  tests 
performed  in  this  work.  For  the  quasi-static  and  dynamic  simulations,  the 
modified  HJC  provided  better  results  for  damaged  material  when  com¬ 
pared  to  the  experimental  tests. 

An  analysis  of  the  effects  of  inertia  and  friction  on  the  dynamic  strength  of 
the  material  during  the  SHPB  test  was  also  performed.  The  analysis  indi¬ 
cated  that  inertia  did  not  increase  the  dynamic  strength  of  the  material 
likely  due  to  the  relatively  low  strain  rates  used.  Confinement  due  to  fric¬ 
tion  between  the  bar-sample  interfaces  did  not  cause  a  significant  increase 
in  the  sample  strength.  However,  the  analysis  showed  that  effective  confin¬ 
ing  pressures  of  10  MPa  was  possible  during  the  SHPB  tests  for  a  constant 
friction  coefficient  of  0.1.  Because  10  MPa  confining  pressures  eliminated 
the  effects  of  microcracking  from  the  quasi-static  strength  of  the  material, 
the  presence  of  pressures  of  this  magnitude  prompts  further  study  of  the 
effects  of  friction  on  microstructural  mechanisms  during  high  rate  tests. 


ERDC/GSL  TR-17-15 


70 


6  Conclusions 

This  report  investigated  the  effects  of  thermally  induced  microcracking  on 
the  mechanical  properties  of  Salem  limestone.  The  investigation  included 
three  parts:  introduction  of  quantifiable  levels  of  microcracks  by  means  of 
thermal  treating,  mechanical  testing  of  limestone  samples  with  varying 
levels  of  microcracks,  and  modification  of  a  commonly  used  numerical 
model  to  incorporate  the  measured  effects.  This  work  demonstrated  that 
the  three  part  approach  is  useful  for  examination  of  the  mechanical  effects 
of  microcracking  on  quasi-brittle  materials.  It  also  demonstrates  that  the 
information  can  be  used  to  improve  the  predictive  capabilities  of  material 
models  incorporating  the  effects  of  microcrack  damage. 

Salem  limestone  was  chosen  because  of  its  homogeneity  and  its  behavior  is 
similar  to  that  of  other  quasi-brittle  materials  such  as  concrete  and  con¬ 
crete  mortars.  Its  homogeneity  simplified  the  microscale  quantification  by 
isolating  microcracking  as  the  only  microscale  mechanism  involved  in  the 
material  characterization.  The  limestone’s  homogeneous  nature  also  en¬ 
hanced  the  determination  of  microcrack  effects  by  allowing  for  consist¬ 
ently  reproducible  results  in  the  mechanical  tests. 

Thermal  treating  was  found  to  induce  quantifiable  levels  of  microcracks  in 
a  consistent,  reproducible  way  in  Salem  limestone.  Sample  preparation  us¬ 
ing  thermal  treating  and  non-destructive  evaluation  were  economical  and 
straightforward.  Using  this  method,  multiple  test  samples  can  be  produced 
simultaneously  in  a  period  of  about  l  day.  Induced  damage  was  controlled 
primarily  by  one  experimental  parameter,  maximum  temperature.  The 
range  of  maximum  temperatures  required  to  induce  various  levels  of  dam¬ 
age  was  largely  making  correlation  with  specific  damage  levels  straight  for¬ 
ward,  and  reproducibility  of  specific  damage  levels  consistent. 

Non-destructive  evaluation  of  treated  samples  using  sonic  wave  velocities 
was  a  quick  and  effective  means  for  assessing  desired  damage  levels.  The 
sonic  wave  testing  takes  only  a  few  minutes  for  each  sample  before  and  af¬ 
ter  the  heat  treatments.  The  changes  in  sonic  wave  velocities,  which  de¬ 
crease  with  increasing  maximum  heat  treatment  temperatures,  indicated 
that  the  induced  microstructural  changes  were  a  function  of  the  maximum 
temperature.  The  wave  velocities  showed  little  variability  for  a  particular 
choice  of  maximum  heat  treatment  temperature  that  demonstrated  the  ef¬ 
fectiveness  of  the  approach  for  inducing  consistent  levels  of  microcracking. 
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XRD,  DSC,  and  TGA  confirmed  that  no  composition  changes  occurred  at 
500  °C  or  below;  therefore,  room  temperature  to  500  °C  became  the  tem¬ 
perature  range  of  interest  for  the  microcrack  characterization  and  me¬ 
chanical  testing. 

CT  scanning,  SEM,  and  OM  were  used  to  directly  observe  microstructural 
changes  caused  by  the  heat  treatments.  The  OM  analysis  was  found  to  be 
the  most  useful  method  for  microcrack  characterization  yielding  both 
qualitative  and  quantitative  data. 

OM  images  showed  an  increase  in  grain  boundary  and  intragranular 
cracking  with  increasing  maximum  temperatures  in  the  heat  treatments. 
Stereological  evaluation  provided  microcrack  densities  for  the  correspond¬ 
ing  heat  treatments  and  indicated  that  microcrack  density  increased  as 
function  of  the  maximum  temperature  in  the  heat  treatments. 

Quasi-static  testing  was  performed  to  characterize  the  mechanical  re¬ 
sponse  of  the  intact  and  damaged  limestone.  The  quasi-static  tests  in¬ 
cluded  uniaxial  compression,  triaxial  compression,  hydrostatic 
compression,  and  uniaxial  strain/ constant  volume  tests.  These  tests 
demonstrated  that  microcracking  lowers  peak  strengths  (by  40  percent  for 
the  high  damage  level)  at  confining  pressures  lower  than  10  MPa.  Mi¬ 
crocracking  did  not  affect  the  limestone’s  strength  at  pressures  greater 
than  10  MPa.  The  initial  shear  modulus  was  reduced  from  approximately 
14,000  MPa  to  4,000  MPa  (for  UC  tests).  This  change  was  observed  for  all 
levels  of  confining  pressure,  but  the  effects  decreased  with  increasing  con¬ 
fining  pressure. 

Unconfined  dynamic  compression  tests  were  also  performed  on  undam¬ 
aged  and  intact  limestone  using  a  modified  SHPB.  These  tests  showed  that 
microcracks  lower  the  stiffness  of  the  material  even  at  higher  strain  rates. 
In  contrast  to  the  quasi-static  testing,  initial  microcrack  density  did  not 
have  an  effect  on  the  dynamic  strength  of  the  limestone.  Peak  strengths  of 
the  limestone  in  the  quasi-static  tests  were  approximately  72  MPa;  peak 
strengths  in  the  dynamic  tests  were  approximately  90  MPa  regardless  of 
the  damage  level.  Though  similar  incident  stress  pulses  were  imparted  to 
the  samples,  higher  strain  rates  were  observed  for  increasing  levels  of 
damage.  This  is  most  likely  due  to  an  initial  lower  stiffness  in  the  samples 
(due  to  cracks  closing),  which  offers  little  resistance  initially  that  allows 
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the  strain  rate  to  increase  at  a  higher  rate.  This  behavior  may  be  responsi¬ 
ble  for  undamaged  and  damage  samples  having  the  same  dynamic 
strengths.  A  limited  number  of  samples  were  available  for  SHPB  tests,  and 
more  dynamic  tests  are  needed  to  confirm  if  these  trends  remain  for  a 
larger  sample  size. 

The  mechanical  responses  of  the  intact  and  damaged  limestone  were  used 
to  modify  the  HJC  model.  Modifications  were  made  to  account  for  shear 
modulus  degradation  and  failure  surface  changes,  which  were  the  main  ef¬ 
fects  of  damage  on  the  material  behavior  observed  in  the  experimental 
data.  The  original  and  modified  HJC  models  were  used  in  a  numerical 
analysis  of  the  quasi-static  and  dynamic  experiments  performed  in  this 
work.  For  the  quasi-static  and  dynamic  loading,  the  modified  HJC  pro¬ 
vided  better  results  for  damaged  material  when  compared  with  the  experi¬ 
mental  tests. 

An  analysis  of  the  effects  of  inertia  and  friction  on  the  dynamic  strength  of 
the  material  during  the  SHPB  test  was  also  performed.  The  analysis  pro¬ 
vided  evidence  that  inertia  and  friction  did  not  increase  the  dynamic 
strength  of  the  material  likely  due  to  the  relatively  low  strain  rates  used. 
However,  the  analysis  showed  that  effective  confining  pressures  of  10  MPa 
was  possible  during  the  SHPB  tests  for  a  constant  friction  coefficient  of 
o.i.  Because  10  MPa  confining  pressures  eliminated  the  effects  of  mi¬ 
crocracking  from  the  quasi-static  strength  of  the  material,  the  presence  of 
pressures  of  this  magnitude  prompts  further  study  of  the  effects  of  friction 
on  microstructural  mechanisms  during  high  rate  tests  especially  if  strain 
rates  greater  than  100  s_1  are  used. 

In  summary,  this  report  demonstrates  a  systematic  approach,  that  is  both 
simple  and  economical,  to  evaluate  microcracking  and  its  effects  on  the 
mechanical  response  of  a  homogeneous,  quasi-brittle  material.  This  ap¬ 
proach  provides  a  starting  point  for  examining  microcracking  in  other 
quasi-brittle  materials  (e.g.,  concrete  mortars,  simple  concrete  compo¬ 
sites)  in  which  microcracking  is  one  of  the  fundamental  components  of  the 
failure  process.  Microcracks  are  caused  by  many  processes  including  me¬ 
chanical  loading,  thermal  loading,  freezing/thawing,  and  chemical  pro¬ 
cesses  (especially  alkali-silica  reactions)  and  understanding  its  effects  on 
these  materials  is  vital  because  quasi-brittle  materials  are  universally  used 
in  military  structures  and  civil  infrastructure.  This  work  also  provides  the 
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basis  for  more  rapid  testing  of  damaged  materials,  which  is  vital  in  evalu¬ 
ating  the  mechanical  response  of  a  structure  subjected  to  microcracking. 

The  mechanical  test  data  provided  by  the  procedures  developed  in  this 
work  is  also  valuable  in  the  development  of  predictive  models.  The  experi¬ 
mental  data  obtained  with  these  procedures  was  used  to  establish  basic 
principles  of  damaged  limestone  behavior  that  were  implemented  in  a 
simple  material  model  to  provide  a  more  realistic  response  to  damage.  Ap¬ 
plying  this  systematic  approach  to  other  quasi-brittle  materials,  such  as 
concrete,  would  provide  much  needed  experimental  data  that  could  be  im¬ 
plemented  into  both  simple  and  complex  material  models.  Improved  dam¬ 
age  effects  implementation  is  valuable  in  the  design  and  analysis 
capabilities  of  structures  composed  of  quasi-brittle  materials. 
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Appendix  A:  Quasi-Static  Testing  Results 

Figure  A.l.  Principal  stress  difference  vs.  axial  strain  for  all  unconfined  compression 

tests. 
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Figure  A.2.  Principal  stress  difference  vs.  volumetric  strain  for  all  unconfined 

compression  tests. 
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Figure  A.3.  Principal  stress  difference  vs.  axial  strain  for  all  triaxial  compression  tests 

with  10  MPa  confining  pressure. 


Figure  A.4.  Principal  stress  difference  vs.  volumetric  strain  for  all  triaxial  compression 
tests  with  10  MPa  confining  pressure. 
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Figure  A.5.  Principal  stress  difference  vs.  axial  strain  for  all  triaxial  compression  tests 

with  20  MPa  confining  pressure. 
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Figure  A.6.  Principal  stress  difference  vs.  volumetric  strain  for  all  triaxial  compression 
tests  with  20  MPa  confining  pressure. 
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Figure  A. 7.  Principal  stress  difference  vs.  axial  strain  for  all  triaxial  compression  tests 

with  100  MPa  confining  pressure. 


Figure  A.8.  Principal  stress  difference  vs.  volumetric  strain  for  all  triaxial  compression 
tests  with  100  MPa  confining  pressure. 
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Figure  A.9.  Principal  stress  difference  vs.  axial  strain  for  all  triaxial  compression  tests 

with  400  MPa  confining  pressure. 


Figure  A.10.  Principal  stress  difference  vs.  volumetric  strain  for  all  triaxial  compression 
tests  with  400  MPa  confining  pressure. 
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Figure  A.ll.  Axial  stress  vs.  axial  strain  for  \JX/CV  tests  with  UX  phase  up  to  100  MPa 

pressure. 


Figure  A.12.  Mean  normal  stress  vs.  volumetric  strain  for  UX/CV  tests  with  UX  phase 

up  to  100  MPa  pressure. 
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Figure  A.13.  Principal  stress  difference  vs.  mean  normal  stress  for  UX/CV  tests  with 

UX  phase  up  to  100  MPa  pressure. 
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Figure  A.14.  Principal  stress  difference  vs.  axial  strain  for  UX/CV  tests  with  UX  phase 

up  to  100  MPa  pressure. 
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Figure  A.15.  Axial  stress  vs.  axial  strain  for  \JX/CV  tests  with  UX  phase  up  to  150  MPa 

pressure. 
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Figure  A.16.  Mean  normal  stress  vs.  volumetric  strain  for  UX/CV  tests  with  UX  phase 

up  to  150  MPa  pressure. 
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Figure  A.17.  Principal  stress  difference  vs.  mean  normal  stress  for  UX/CV  tests  with 

UX  phase  up  to  150  MPa  pressure. 
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Figure  A.18.  Principal  stress  difference  vs.  axial  strain  for  UX/CV  tests  with  UX  phase 

up  to  150  MPa  pressure. 
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Figure  A.19.  Axial  stress  vs.  axial  strain  for  UX/CV  tests  with  UX  phase  up  to  200  MPa 

pressure. 
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Figure  A.20.  Mean  normal  stress  vs.  volumetric  strain  for  UX/CV  tests  with  UX  phase 

up  to  200  MPa  pressure. 
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Figure  A.21.  Principal  stress  difference  vs.  mean  normal  stress  for  UX/CV  tests  with 

UX  phase  up  to  200  MPa  pressure. 
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Figure  A.22.  Principal  stress  difference  vs.  axial  strain  for  UX/CV  tests  with  UX  phase 

up  to  200  MPa  pressure. 
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Figure  A.23.  Mean  normal  stress  vs.  volumetric  strain  for  HC  tests. 

350 


300 

(0 

CL 

S.  250 

a> 

o 

c 

Si 

£  200 

b 

a) 

</) 

S  150 
c/) 

i  ioo 

‘C 

Q- 


50 


0 

0  0.02  0.04  0.06  0.08  0.1  0.12  0.14  0.16 

Axial  Strain 


ERDC/GSL  TR-17-15 


90 


Appendix  B:  Dynamic  Testing  Results 


Figure  B.l.  Interface  stresses  and  strain  rate  from  modified  SHPB  test  on  undamaged 

sample  (sample  12F). 
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Figure  B.2.  Interface  stresses  and  strain  rate  from  modified  SHPB  test  on  undamaged 

sample  (sample  15D). 
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Figure  B.3.  Interface  stresses  and  strain  rate  from  modified  SHPB  test  on  low  damage 

sample  (sample  17D). 
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Figure  B.4.  Interface  stresses  and  strain  rate  from  modified  SHPB  test  on  low  damage 

sample  (sample  17B). 


Figure  B.5.  Interface  stresses  and  strain  rate  from  modified  SHPB  test  on  low  damage 

sample  (sample  22A). 
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Figure  B.6.  Interface  stresses  and  strain  rate  from  modified  SHPB  test  on  high 
damage  sample  (sample  17C). 
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Figure  B.8.  Axial  stress  vs.  axial  strain  for  all  SHPB  tests. 
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Appendix  C:  Mechanical  Test  Sample 
Properties 


Table  C.l.  Mechanical  test  sample  properties. 


Test 

Number 

Type  of  Test 

ID 

Height,  in 

Dia,  in 

wt,g 

Rad  Stress,  MPa 

Wet  Density, 
Mg/m3 

llslaOl 

UC 

SL-01-J-01 

4.382 

2.000 

534.73 

2.370 

llsla02 

UC 

SL-01-J-02 

4.444 

2.000 

542.30 

2.370 

llsla03 

UC 

SL-01-K-01 

4.360 

2.001 

532.83 

2.371 

llsla04 

UC 

SL-01-G-01 

4.429 

2.013 

540.77 

2.341 

1 1  sla05 

UC 

SL-01-G-02 

4.463 

2.012 

543.77 

2.339 

1 1  sla06 

UC 

SL-0 1-1-02 

4.439 

2.013 

541.66 

2.340 

llsla07 

UC 

SL-01-H-01 

4.424 

2.017 

532.39 

2.298 

1 1  sla08 

UC 

SL-01-H-02 

4.443 

2.018 

534.09 

2.294 

llsla09 

UC 

SL-01-I-01 

4.391 

2.021 

535.55 

2.320 

llslbOl 

TXC/20 

SL-01-T-01 

4.344 

2.008 

519.60 

20 

2.305 

1 1  slb02 

TXC/20 

SL-0 1-0-01 

4.416 

1.999 

535.46 

20 

2.358 

1 1  slb03 

TXC/20 

SL-01-Q-01 

4.439 

2.000 

542.51 

20 

2.374 

1 1  slb04 

TXC/20 

SL-01-U-01 

4.413 

2.008 

536.91 

20 

2.344 

1 1  slb05 

TXC/100 

SL-0 1 -P-01 

4.402 

1.999 

543.22 

100 

2.399 

1 1  slb06 

TXC/100 

SL-01-R-02 

4.370 

2.000 

527.09 

100 

2.343 

1 1  slb07 

TXC/100 

SL-01-U-02 

4.412 

2.008 

534.62 

100 

2.335 

llslb08 

TXC/100 

SL-01-T-02 

4.438 

2.008 

530.45 

100 

2.303 

llslb09 

TXC/20 

SL-01-Y-01 

4.418 

2.002 

523.22 

20 

2.296 

llslblO 

TXC/20 

SL-01-Y-02 

4.384 

2.002 

528.93 

20 

2.339 

11  slbll 

TXC/100 

SL-01-X-02 

4.344 

2.003 

539.48 

100 

2.405 

1 1  sib  1 2 

TXC/100 

SL-01-Z-01 

4.338 

2.001 

543.27 

100 

2.430 

llslbl3 

HC 

SL003 

4.374 

2.004 

526.80 

400 

2.330 

1 1  sib  1 4 

HC 

SL014 

4.410 

2.012 

533.19 

400 

2.321 

1 1  sib  1 5 

HC 

SL026 

4.354 

2.008 

528.19 

400 

2.338 

1 1  sib  1 6 

TXC/10 

SL006 

4.396 

2.004 

531.65 

10 

2.340 

1 1  sib  1 8 

TXC/10 

SL015 

4.374 

2.017 

535.84 

10 

2.340 

1 1  sib  1 9 

TXC/10 

SL016 

4.409 

2.019 

536.52 

10 

2.319 

1 1  slb20 

TXC/10 

SL029 

4.455 

2.009 

540.93 

10 

2.337 

1 1  slb2 1 

TXC/400 

SL001 

4.404 

2.003 

529.43 

10 

2.328 

llslb22 

TXC/10 

SL030 

4.419 

2.008 

537.46 

10 

2.344 

llslb23 

TXC/400 

SL004 

4.420 

2.003 

532.23 

400 

2.332 

1 1  slb24 

TXC/400 

SL022 

4.467 

2.017 

538.63 

400 

2.303 

llslb25 

TXC/400 

SL023 

4.438 

2.016 

532.60 

400 

2.294 

llslb26 

TXC/400 

SL024 

4.460 

2.009 

538.25 

400 

2.323 

1 1  slb27 

TXC/400 

SL025 

4.455 

2.009 

539.22 

400 

2.330 

1 1  slb28 

UX/CV/100 

SL008 

4.269 

2.003 

516.30 

100 

2.342 

1 1  slb29 

UX/CV/100 

SL017 

4.411 

2.018 

537.05 

100 

2.323 

llslb30 

UX/CV/100 

SL027 

4.464 

2.009 

543.41 

100 

2.343 

1  lslb3 1 

UX/CV/150 

SL009 

4.442 

2.004 

536.79 

150 

2.338 

1 1  slb32 

UX/CV/150 

SL018 

4.469 

2.018 

544.79 

150 

2.326 

llslb33 

UX/CV/150 

SL028 

4.404 

2.011 

533.39 

150 

2.327 

1 1  slb34 

UX/CV/200 

SL011 

4.355 

2.003 

526.08 

200 

2.339 

llslb35 

UX/CV/200 

SL020 

4.443 

2.018 

540.61 

200 

2.322 

llslb36 

UX/CV/200 

SL032 

4.276 

2.019 

517.91 

200 

2.309 

llslb37 

TXC/400 

SL-01-X-01 

4.358 

2.004 

525.23 

400 

2.332 
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Table  C.2.  Additional  mechanical  test  sample  properties. 


Test  Number 

Specimen  ID 

Posttest 
Water 
Content,  % 

Dry 

Density, 

Mg/m3 

Porosity,  % 

Degree  of 
Saturation,  % 

Volume  of 
Air,  % 

Volume  of 
Water,  % 

Volume  of 
Solids,  % 

llslaOl 

SL-01-J-01 

0.0 

2.369 

13.75 

0.86 

13.64 

0.12 

86.25 

1 1  sla02 

SL-01-J-02 

0.0 

2.369 

13.76 

1.03 

13.62 

0.14 

86.24 

1 1  sla03 

SL-01-K-01 

0.06 

2.370 

13.72 

1.04 

13.58 

0.14 

86.28 

1 1  sla04 

SL-01-G-01 

0.04 

2.340 

14.81 

0.63 

14.71 

0.09 

85.19 

1 1  sla05 

SL-01-G-02 

0.04 

2.338 

14.90 

0.63 

14.81 

0.09 

85.10 

1  lsla06 

SL-0 1-1-02 

0.04 

2.339 

14.86 

0.63 

14.77 

0.09 

85.14 

1 1  sla07 

SL-01-H-01 

0.03 

2.298 

16.36 

0.42 

16.29 

0.07 

83.64 

1 1  sla08 

SL-01-H-02 

0.03 

2.293 

16.53 

0.42 

16.46 

0.07 

83.47 

1 1  sla09 

SL-01-I-01 

0.04 

2.319 

15.57 

0.60 

15.48 

0.09 

84.43 

llslbOl 

SL-01-T-01 

0.05 

2.304 

16.13 

0.71 

16.02 

0.12 

83.87 

1 1  slb02 

SL-0 1-0-01 

0.1 

2.355 

14.26 

1.65 

14.02 

0.24 

85.74 

1 1  slb03 

SL-01-Q-01 

0.1 

2.372 

13.67 

1.74 

13.43 

0.24 

86.33 

1 1  slb04 

SL-01-U-01 

0.0 

2.344 

14.69 

0.64 

14.59 

0.09 

85.31 

1 1  slb05 

SL-0 1 -P-01 

0.18 

2.395 

12.81 

3.37 

12.38 

0.43 

87.19 

1 1  slb06 

SL-01-R-02 

0.18 

2.339 

14.86 

2.83 

14.44 

0.42 

85.14 

1 1  slb07 

SL-01-U-02 

0.0 

2.333 

15.06 

1.08 

14.89 

0.16 

84.94 

1  lslb08 

SL-01-T-02 

0.07 

2.302 

16.21 

0.99 

16.05 

0.16 

83.79 

1 1  slb09 

SL-01-Y-01 

0.1 

2.293 

16.52 

1.67 

16.25 

0.28 

83.48 

llslblO 

SL-01-Y-02 

0.1 

2.336 

14.97 

2.03 

14.66 

0.30 

85.03 

11  slbll 

SL-01-X-02 

0.0 

2.403 

12.51 

1.35 

12.34 

0.17 

87.49 

1 1  sib  1 2 

SL-01-Z-01 

0.06 

2.429 

11.59 

1.26 

11.44 

0.15 

88.41 

1 1  sib  1 3 

SL003 

0.12 

2.327 

15.28 

1.83 

15.00 

0.28 

84.72 

1 1  sib  1 4 

SL014 

0.07 

2.319 

15.58 

1.04 

15.42 

0.16 

84.42 

1 1  sib  1 5 

SL026 

0.18 

2.333 

15.05 

2.79 

14.63 

0.42 

84.95 

1 1  sib  1 6 

SL006 

0.05 

2.339 

14.87 

0.79 

14.75 

0.12 

85.13 

1 1  sib  1 8 

SL015 

0.03 

2.339 

14.85 

0.47 

14.78 

0.07 

85.15 

1 1  sib  1 9 

SL016 

0.04 

2.319 

15.60 

0.59 

15.51 

0.09 

84.40 

1 1  slb20 

SL029 

0.07 

2.336 

14.97 

1.09 

14.80 

0.16 

85.03 

1 1  slb2 1 

SL001 

0.05 

2.327 

15.29 

0.76 

15.17 

0.12 

84.71 

1 1  slb22 

SL030 

0.07 

2.342 

14.74 

1.11 

14.58 

0.16 

85.26 

1 1  slb23 

SL004 

0.17 

2.328 

15.25 

2.59 

14.86 

0.40 

84.75 

1 1  slb24 

SL022 

0.08 

2.301 

16.23 

1.13 

16.05 

0.18 

83.77 

1  lslb25 

SL023 

0.09 

2.292 

16.56 

1.25 

16.35 

0.21 

83.44 

1 1  slb26 

SL024 

0.20 

2.319 

15.59 

2.97 

15.13 

0.46 

84.41 

1 1  slb27 

SL025 

0.20 

2.325 

15.35 

3.03 

14.88 

0.47 

84.65 

1  lslb28 

SL008 

0.12 

2.339 

14.84 

1.89 

14.56 

0.28 

85.16 

1 1  slb29 

SL017 

0.06 

2.322 

15.49 

0.90 

15.35 

0.14 

84.51 

1  lslb30 

SL027 

0.15 

2.340 

14.82 

2.37 

14.47 

0.35 

85.18 

1 1  slb3 1 

SL009 

0.12 

2.335 

14.99 

1.87 

14.71 

0.28 

85.01 

1  lslb32 

SL018 

0.06 

2.324 

15.38 

0.91 

15.24 

0.14 

84.62 

1  lslb33 

SL028 

0.18 

2.323 

15.44 

2.71 

15.03 

0.42 

84.56 

1  lslb34 

SL011 

0.12 

2.337 

14.94 

1.88 

14.66 

0.28 

85.06 

1  lslb35 

SL020 

0.07 

2.320 

15.55 

1.04 

15.39 

0.16 

84.45 

1  lslb36 

SL032 

0.19 

2.304 

16.12 

2.72 

15.68 

0.44 

83.88 

1  lslb37 

SL-01-X-01 

0.1 

2.328 

15.24 

2.29 

14.90 

0.35 

84.76 
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Appendix  D:  SHPB  Simulation  Input  File 


$ 

$ 

$ 

UNITS 

pressure : 

dyne/cmA2 

$ 

length 

cm 

$ 

time 

second 

$ 

mass 

gram 

$ 

force 

dyne 

$ 

Temp. 

degK 

$ 

Velocity 

cm/ sec 

$ 

Density 

g/cc 

$ 

Force 

dyne 

$ 

Energy 

erg 

$ 

Power 

erg/ sec 

$ 

10000000 

$  1  MPa 

$ 

100000 

$  1  newton 

$ 

le+12 

$  1  Mbar 

$ 

0.001 

$  1  ms 

$ 

100 

$  1  meter/sec 

$ 

1 

$  1  g/cc 

$ 

2 

$  cube  root  of  8 

begin  sierra 

begin  definition  for  function  12 Fstress_wave 
type  is  piecewise  linear 
begin  values 


0.0 

0. 

0 

o 

CO 

s-6 

3. 

6e7 

13.' 

Oe-6 

12 

.  9e7 

20. 

5e-6 

27 

.  9e7 

75.' 

Oe-6 

90 

.  0e7 

100 

.  5e- 

6 

11 

4 . 8e7 

108 

.  5e- 

6 

11 

8 .  Ie7 

115 

.  Oe- 

6 

11 

8 . 8e7 

127 

.  Oe- 

6 

11 

5 . 7e7 

135 

.  Oe- 

6 

11 

0 . 4e7 

166 

.0e- 

6 

80 

.  4e7 

174 

.  5e- 

6 

67 

.  2e7 

190 

.0e- 

6 

33 

.  3e7 

200 

.  Oe- 

6 

18 

.  6e7 

210 

.  Oe- 

6 

9. 

3e7 

215 

.  5e- 

6 

6. 

0e7 

227 

.  5e- 

6 

1. 

8e7 

232 

.  Oe- 

6 

0. 

0 

end  values 

end 

begin  definition  for  function  15Dstress_wave 
type  is  piecewise  linear 
begin  values 

0.0  0.0 
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7.0e-6  1.5e7 

13 . 5e-6  9 . 3e7 

20 . 5e-6  27 . 9e7 

96 . 5e-6  121. 4e7 

109 . 5e-6  128. 6e7 
117 . 5e-6  128. 6e7 
12  6 . Oe-6  125. 6e7 
140 . Oe-6  116. 0e7 
14  9 . 5e-6  107. 6e7 
17 2 . 5e-6  80 . 3e7 
183. Oe-6  54 . 6e7 
192. Oe-6  36 . 3e7 
201 . 5e-6  22 . 5e7 
209. 5e-6  15.0e7 
223. Oe-6  8 . Ie7 
238 . 5e-6  5.1e7 
253 . 5e-6  3.0e7 
266. Oe-6  0.0 
end  values 

end 

begin  definition  for  function  damaged_stress_wave 
type  is  piecewise  linear 
begin  values 

0.0  0.0 
15. Oe-6  27 . 0e7 
30. Oe-6  44 . 4e7 
81 . 5e-6  104. 4e7 

100. Oe-6  122. Ie7 
112. Oe-6  124. 4e7 
118. Oe-6  123. 0e7 
133. Oe-6  112. 5e7 
167 . 5e-6  7 1 . 7e7 
195. Oe-6  16 . 8e7 
211. Oe-6  3 . 0e7 
220 . 5e-6  0.0 
end  values 

end 

begin  property  specification  for  material  steel 
density  =  8.0 

begin  parameters  for  model  elastic 
youngs  modulus  =  2.00e+12 
poissons  ratio  =  0.35 
end  parameters  for  model  elastic 

end  property  specification  for  material  steel 

begin  property  specification  for  material  salem_lime 
density  =  2.3 
bulk  modulus  =  15000e7 
shear  modulus  =  10000e7 


$13000e7 
$  98  60e7 
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begin  parameters  for  model  h j cconcrete_mmm 
init  density  =  2.3 
init  shear  modulus  =  10000e7 
yield  stress  a  =  0.50 

press  hard  coef  b  =  1.286  $1,373  $1,428  no  fric 

$1,286  no  rate 

press  hard  exp  n  =  0.8362  $0.8273  $0.8223  no  fric 

$0.8362  no  rate 

strain  rate  coef  c  =  0.0  $0.0038  $6.0e-3 

comp  stren  fc  =  72.0e7 

max  stress  =  13.0 

pcrush  =  195. 0e7 

ucrush  =  0 . 013 

press  coef  kl  =  44300. 0e7 

press  coef  k2  =  0.0 

press  coef  k3  =  0.0 

plocki  =  665. 0e7 

ulock  =  0.134 

max  tens  press  t  =  6 . 0e7 

damage  coef  dl  =  0.05 

damage  exp  d2  =  1.0 

min  fail  strain  =  0.011 

art  vis  cl  =  0.2 

art  vis  cq  =  4.0 

freeOl  =  1.0 

f ree02  =  0.65 

free03  =  1.0 

free04  =3.5 

f ree05  =  8.0 

free06  =  0.0 

free07  =0.0 

free08  =0.0 

free09  =  1.0 

freelO  =  2000. 0e7 

end  parameters  for  model  h j cconcrete_mmm 

end  property  specification  for  material  salem_lime 

begin  finite  element  model  SHPB 

database  name  =  CAVS_SHPB_cm_f ine2 . g 
$  database  type  =  exodusii' 

begin  parameters  for  block  block_l 
material  steel 

solid  mechanics  use  model  elastic 
end  parameters  for  block  block_l 

begin  parameters  for  block  block_2 
material  steel 

solid  mechanics  use  model  elastic 
end  parameters  for  block  block_2 
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begin  parameters  for  block  block_3 
material  salem_lime 

solid  mechanics  use  model  h j cconcrete_mmm 
end  parameters  for  block  block_3 
end  finite  element  model  SHPB 

begin  presto  procedure  presto_procedure 

begin  time  control 

begin  time  stepping  block  stress_input 
start  time  =  0.0 

begin  parameters  for  presto  region  presto_region 
$  step  interval  =  25 

$  time  step  scale  factor  =  0.10 

end  parameters  for  presto  region  presto_region 
end  time  stepping  block  stress_input 
termination  time  =  1000e-6 
end  time  control 

begin  presto  region  presto_region 
use  finite  element  model  SHPB 

$  begin  results  output  results 

$  database  name  =  hoppyBar_CAVS . e 

$  At  Time  0.0  Increment  =  l.e-08 

$  nodal  Variables  =  displacement 

$  nodal  Variables  =  velocity 

$  element  Variables  =  log_strain 

$  element  Variables  =  stress 

$  element  Variables  =  von_mises 

$  global  Variables  =  timestep 

$  global  Variables  =  time 

$  end  results  output  results 

begin  results  output  sample_results 
database  name  =  hjc_shpb.e 
at  time  0.0  increment  =  5.e-06 
nodeset  variables  =  displacement 
nodeset  variables  =  velocity 
element  variables  =  failure_flag 
element  variables  =  log_strain 
element  variables  =  stress 
element  variables  =  von_mises 
element  variables  =  eqps 
element  variables  =  damage 
element  variables  =  pressure 
global  variables  =  timestep 
global  variables  =  time 
start  time  =  400.e-06 
termination  time  =  600.e-06 
end  results  output  sample_result s 
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begin  heartbeat  output  stress_strain 
stream  name  =  stress_strain.txt 
start  time  =  0.0 
at  time  0.0  increment  =  5.0e-07 
element  log_strain (xx)  nearest  location  114, 
0.452548,  0.452548  as  log_strain_01 

element  log_strain (xx)  nearest  location  287, 
0.452548,  0.452548  as  log_strain_02 
format  =  default 
precision  =  10 
labels  =  off 
legend  =  on 
timestamp  format  "" 
end  heartbeat  output  stress_strain 


begin  heartbeat  output  stress_strain 

stream  name  =  element_stressStrain.txt 

start  time  =  0.0 

at  time  0.0  increment  =  5.0e-07 

element  log_strain (xx)  at  element  85335  as 


ax  strian  cent 


element  stress (xx)  at  element  85355  as  ax_str_cent 
element  pressure  at  element  85335  as  press_cent 
element  damage  at  element  85335  as  dam_cent 
element  eqps  at  element  85335  as  eqps_cent 
element  log_strain (xx)  at  element  85361  as 


ax_strain_edge 


element  stress  (xx)  at  element  85361  as 


ax_stress_edge 


element  pressure  at  element  85361  as  press_edge 
element  damage  at  element  85361  as  dam_edge 
element  eqps  at  element  85361  as  eqps_edge 
element  log_strain (xx)  at  element  85347  as 


$ 

$ 

$ 

$ 

ax_strain_mid 

$  element  stress (xx)  at  element  85347  as  ax_stress_mid 

$  element  pressure  at  element  85347  as  press_mid 

$  element  damage  at  element  85347  as  dam_mid 

$  element  eqps  at  element  85347  as  eqps_mid 

$  precision  =  10 

$  format  =  default 

$  labels  =  off 

$  legend  =  on 

$  timestamp  format  " " 

$  end  heartbeat  output  stress_strain 

$$$$$$$$$$$$$$$ 

$$$  BCs  $$$ 

$$$$$$$$$$$$$$$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$  Axisymmetric  BC  for  Bars 

$$$$$$$$$$$$$$$$$$$$$$$$$$ 

begin  fixed  displacement 

active  periods  =  stress_input 
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node  set  =  nodelist_2 
node  set  =  nodelist_6 
node  set  =  nodelist_10 
component  =  y 

end 

begin  fixed  displacement 
node  set  =  nodelist_4 
node  set  =  nodelist_8 
node  set  =  nodelist_12 
component  =  z 

end 

$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$  Stress  Wave  Input 
$$$$$$$$$$$$$$$$$$$$$$$$$$ 
begin  pressure 

surface  =  surface_l 

$  function  =  damaged_stress_wave 

function  =  15Dstress_wave 
scale  factor  =  1 
end  pressure 

$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$  Element  Death 

$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$  begin  element  death 

$  block  =  block_3 

$$  criterion  is  element  value  of  log_strain (yy)  >= 

0.020 

$$  criterion  is  element  value  of  log_strain ( z z )  >= 

0.020 

$$$$  criterion  is  element  value  of 

ef f ective_log_strain  >=  0.020 

$  criterion  is  element  value  of  eqps  >=  0.018 

$$  death  method  =  disconnect  element 

$  end  element  death 

$ 

$  begin  element  death 

$  block  =  block_3 

$  death  on  inversion  =  on 

$  criterion  is  element  value  of  nodal_j acobian_ratio 

<-0.0 

$  end  element  death 

$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$  Contact 

$$$$$$$$$$$$$$$$$$$$$$$$$$ 

begin  contact  definition  incident 

$  contact  surface  incident_bar  contains  block_l 

$  contact  surface  transmitted_bar  contains  block_2 

contact  surface  sample  contains  block_3 

begin  contact  surface  incident_bar 
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surface  =  surface_3 
end  contact  surface  incident_bar 

$ 

$  begin  contact  surface  sample_01 

$  surface  =  surface_5 

$  end  contact  surface  sample_01 

$ 

$  begin  contact  surface  sample_02 

$  surface  =  surface_7 

$  end  contact  surface  sample_02 

$ 

begin  contact  surface  transmitted_bar 
surface  =  surface_9 
end  contact  surface  transmitted_bar 

begin  constant  friction  model  friction 
friction  coefficient  =0.5 
end  constant  friction  model  friction 

begin  interaction  interaction_l 
surfaces  =  incident_bar  sample 
friction  model  =  friction 
end  interaction  interaction__l 

begin  interaction  interaction_2 

surfaces  =  transmitted_bar  sample 
friction  model  =  friction 
end  interaction  interaction_2 

end  contact  definition  incident 
end  presto  region  presto_region 
end  presto  procedure  presto_procedure 


end  sierra 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for 
reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA 
22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display 
a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

June  2017  Final 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

Effects  of  Thermally  Induced  Microcracking  on  the  Quasi-Static  and  Dynamic  Response 
of  Salem  Limestone 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

P2  331691 

Zachary  Kyle  Crosby 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

Geotechnical  and  Structures  Laboratory 

U.S.  Army  Engineer  Research  and  Development  Center 

3909  Halls  Terry  Road 

Vicksburg,  MS  39180-6199 

ERDC/GSL  TR-17-15 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

U.S.  Army  Engineer  Research  and  Development  Center 

Vicksburg,  MS 

USACE 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

The  effects  of  microcracking  on  the  mechanical  properties  of  Salem  limestone  were  investigated  in  three  phases:  introduction  of  quanti¬ 
fiable  levels  of  microcracks  by  thermal  treating,  mechanical  testing  of  limestone  samples  with  varying  levels  of  microcracks,  and  modi¬ 
fication  of  a  numerical  model  to  incorporate  the  measured  effects.  Computed  tomography  scanning,  scanning  electron  microscopy,  and 
optical  microscopy  (OM)  were  used  to  observe  micro  structural  changes  caused  by  the  heat  treatments.  Mechanical  testing  was  per¬ 
formed  to  characterize  the  mechanical  response  of  the  intact  and  damaged  limestone.  Quasi-static  tests  included  uniaxial  compression, 
triaxial  compression,  hydrostatic  compression,  and  uniaxial  strain  /  constant  volume  tests.  Microcracking  did  not  affect  the  limestone’s 
strength  at  pressures  greater  than  10  MPa.  Dynamic  tests  were  performed  using  a  modified  split  Hopkinson  pressure  bar. 

The  results  of  the  mechanical  tests  were  used  to  modify  the  HJC  model.  Modifications  were  made  to  account  for  shear  modulus  degra¬ 
dation  and  failure  surface  changes.  The  original  and  modified  HJC  models  were  used  in  a  numerical  analysis  of  the  mechanical  tests 
performed  in  this  work.  The  modified  HJC  provided  better  results  for  damaged  material  when  compared  with  the  quasi- static  and  dy¬ 
namic  experiments.  This  work  demonstrated  that  this  approach  is  useful  for  examination  of  the  effects  of  microcracking  on  quasi-brittle 
materials  and  can  be  used  to  improve  the  predictive  capabilities  of  material  models. 

15.  SUBJECT  TERMS  (see  reverse) 


16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION 

OF  ABSTRACT 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

S  AR 

Unclassified  Unclassified  Unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


18.  NUMBER  19a.  NAME  OF  RESPONSIBLE 
OF  PAGES  PERSON 

19b.  TELEPHONE  NUMBER  (include 
i  i  a  area  code) 


15.  SUBJECT  TERMS  (concluded) 

Quasi-static  testing 

Dynamic  testing 

Limestone  -  Cracking 

Split  Hopkinson  Pressure  Bar  (SHPB) 

Damage  model 

Microcracking 

Thermal  effects 

Materials  -  Mechanical  properties 

Microstructure 

Fracture  mechanics 

Brittleness 

Numerical  analysis 


